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L  ABSTRACT 


The  Center  for  Theoretical  Geo/Cosmo  Plasma  Physics  was  established  by  the  Air 
Force  Office  of  Scientific  Research  in  1986  through  a  DoD  University  Research  Initiative 
(URI)  Grant  via  keen  national  competition.  The  goal  of  the  center  since  its  inception  has 
been  to  develop  and  maintain  a  program  of  excellence  in  interdisciplinary  space  plasma 
research  involving  the  mutual  interactions  of  collaborating  members  of  a  select  group  of 
space  scientists,  plasma  physicists,  mathematicians  and  numerical  analysts.  During  the 
past  several  years,  under  the  grant  title,  "Frontier  Geoplasma  Research",  members  of  the 
center  have  made  seminal  contributions  to  a  number  of  definitive  research  findings 
related  to  the  phenomena  of  intermittent  plasma  turbulence,  forced  and/or  self-organized 
criticality,  global  acceleration  of  the  solar  wind  and  polar  wind,  sporadic  localized 
reconnections  in  the  magnetotail  and  in  the  auroral  zone,  charged  particle  energization 
through  wave-particle  interactions,  the  black  auroral  curls,  multi-scale  evolutions, 
magnetosphere/ionosphere  coupling,  and  the  theory  of  complexity  in  space  plasmas. 
Some  of  the  results  of  these  research  activities  have  already  found  practical  applications 
toward  the  missions  of  the  United  States  Air  Force,  primarily  through  the  collaborating 
efforts  between  the  center  members  and  members  of  the  research  group  headed  by  Dr. 
J.R.  Jasperse  at  the  Air  Force  Research  Laboratory. 

The  MIT  Center  is  now  recognized  by  their  peers  worldwide  as  a  star  example  of 
a  successful  enterprise  between  a  prestigious  educational  institution  and  an  established 
governmental  research  laboratory,  the  Air  Force  Research  Laboratory  [AFRL].  In 
addition  to  our  collaborating  research  efforts,  our  center  and  AFRL  have  co-sponsored  a 
symposium  series  on  the  "Physics  of  Space  Plasmas"  and  a  workshop  series  on 
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"Theoretical  Geoplasma  Physics",  now  under  the  combined  name,  the  "Cambridge 
Symposia/Workshops  on  Space  Plasma  Physics".  The  center  also  periodically  co¬ 
sponsors  other  scientific  workshops  and  symposia  with  other  institutions  on  various 
important  topics  of  current  space  plasma  research.  Proceedings  of  these  symposia  and 
workshops  have  become  important  research  resources  for  these  current  topics.  For 
example,  one  such  recent  volume  was  devoted  to  the  memory  of  Professor 
Subrahmanyan  Chandrasekhar,  a  Nobel  laureate  and  a  giant  in  the  scientific  fields  of 
plasma  physics,  astrophysics  and  space  sciences.  Other  noteworthy  contributions  include 
the  volume  of  a  collection  of  outstanding  contributions  in  the  modem  physics  of 
complexity  published  as  a  special  issue  entitled  "Forced  and/or  Self-Organiced  Criticality 
[FSOC]  in  Space  Plasmas"  by  the  Journal  of  Atmospheric  and  Solar  Terrestrial  Physics 
[JASTP],  and  a  special  topics  section  on  Transports  Processes  in  Geospace  in  the  journal 
of  Plasma  Physics. 

The  center  has  participated  actively  with  the  graduate  school  of  MIT,  in  a 
Minority  Summer  Research  Program,  to  effect  changes  in  the  realm  of  graduate  scientific 
training  for  minorities.  It  is  also  an  active  participant  in  the  Research  Scientist  Institute 
(RSI)  program  with  the  Washington-based  Center  of  Excellence  in  Education  (CEE)  in 
sponsoring  the  preparation  of  genius-level  high  school  students  recruited  worldwide  in 
preparation  for  them  to  achieve  the  ultimate  of  their  capabilities.  CEE  was  founded  by 
Admiral  Rickover  in  1983  and  its  honorary  trustees  include  such  luminaries  as  former 
president  Jimmy  Carter,  Senator  Joseph  Leiberman,  and  Secretary  Colin  Powell.  In 
addition,  the  MIT  center  has  maintained  an  active  visiting  scientists  program,  which  has 
kept  our  research  program  vibrant  and  up-to-date. 
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During  the  period,  the  members  of  the  center  have  published  24  scientific  papers, 
books  and  proceedings,  and  have  been  invited  to  deliver  27  invited  and  review  lectures  at 
various  international  conferences  and  renowned  educational  and  research  institutions. 
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II.  INTRODUCTION 


For  the  United  States  Air  Force  to  enjoy  its  continuing  success  in  meeting  the 
demands  of  its  mission  heading  into  the  21st  century,  it  must  be  prepared  to  operate  in  the 
"New  Frontiers"  of  space,  particularly  its  ever-changing  turbulent  geoplasma 
environment.  Recognizing  this  importance,  a  Center  of  Excellence  in  Theoretical 
Geo/Cosmo  Plasma  Physics  was  established  (through  keen  national  competition)  at  the 
Massachusetts  Institute  of  Technology  under  the  sponsorship  of  the  Air  Force  Office  of 
Scientific  Research. 

Benefiting  from  strong  interactions  with  the  members  of  the  Air  Force  Research 
Laboratory  at  Hanscom,  MA,  our  Center  has  made  substantive  progress  toward  the  goals 
set  forth  since  its  inception.  The  close  proximity  between  the  two  institutions  has 
allowed  direct  personnel  exchanges  as  well  as  direct  and  meaningful  collaborations  in 
numerous  relevant  research  programs  in  geoplasma  physics.  During  the  past  few  years, 
under  the  grant  name,  "Frontier  Geoplasma  Research",  we  have  succeeded  in  jointly 
developing  the  fundamental  understanding  of  several  multiscale  space  plasma 
phenomena  of  Cutting-edge  importance.  These  include  the  study  of  the  origin  of  broad¬ 
band  low  frequency  turbulence  in  the  auroral  zone  and  the  associated  wave-particle 
interactions  such  as  the  anomalous  acceleration  of  ionospheric  oxygen  ions  into 
magnetospheric  energies,  the  global  acceleration  and  evolution  of  fast  solar  streams  into 
the  distant  helisphere,  the  generation  of  intense  divergent  electric  fields  and  the 
phenomena  of  black  auroral  curls,  the  self-consistent  kinetic  evolution  and  the  dynamic 
understanding  of  the  day-night  asymmetry  of  photoelectron-driven  polar  winds,  the  onset 
of  coherent  structures  and  associated  double  layers  in  the  acceleration  region  of  the 
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ionosphere,  the  study  of  intermittent  turbulence  and  sporadic  localized  reconnections  in 
the  Earth's  magnetotail  and  the  solar  wind,  the  formation  of  ion  conics  and  counter¬ 
streaming  electrons  along  auroral  field  lines,  and  the  general  phenomena  of  complexity 
involving  the  coupling  of  the  ionosphere,  magnetosphere,  and  the  solar  wind.  Together, 
we  have  published  (during  the  past  several  years)  24  joint  technical  papers  and  books, 
and  delivered  27  invited  and  review  lectures  at  various  international  conferences  and 
renowned  institutions. 

A  number  of  the  afore-mentioned  research  findings  have  already  found  practical 
applications  by  our  colleagues  at  the  Air  Force  Research  Laboratory.  Examples  of  such 
applications  are:  the  prediction  of  charged-particle  precipitation  patterns  and  deposition 
profiles  in  the  diffuse-auroral  zone  of  the  ionosphere,  the  prediction  of  solar  EUV  and  X- 
ray  fluxes  based  on  the  ionospheric  photoelectron  measurements  and  transport 
calculations,  the  calculation  of  ionospheric  electron  density  profiles  in  the  equatorial, 
mid-,  and  high-latitude  portions  of  the  globe  and  associated  local,  nonlocal,  and  nonlinear 
instabilities,  and  a  fully  self-consistent  analytical  model  that  is  capable  of  accurately 
predicting  the  various  features  of  the  turbulent  geoplasma  environment  (particularly  in 
the  downward  auroral  current  regions)  relevant  to  important  Air  Force  operations. 

Our  Center  and  the  Air  Force  Research  Laboratory  (led  by  Drs.  Tom  Chang  at 
MIT  and  J.R.  Jasperse  at  AFRL)  jointly  inaugurated  a  series  of  symposia  on  the  "Physics 
of  Space  Plasmas".  The  original  purpose  of  these  symposia  is  to  provide  an  annual  get- 
together  for  space  scientists  in  the  Boston-New  England  area.  Its  coverage  has  since 
attracted  contributions  from  scientists  worldwide.  Each  symposium  included  the 
presentation  of  an  "Alfven  lecture"  established  in  honor  of  the  Nobel  Laureate,  Professor 
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Hannes  Alfven  of  the  Swedish  Royal  Institute  of  Technology.  Professor  Alfven  has 
strongly  supported  the  activities  of  the  MIT  center.  He  was  a  founding  honorary  member 
of  the  center  and  visits  the  center  often  during  its  formative  years.  Alfven  Lecturers  have 
included  such  luminaries  as  Professor  Oscar  Buneman  of  the  Stanford  University, 
Professor  Jim  Dungey  of  the  Imperial  College,  London,  Professor  Eugene  Parker  of  the 
University  of  Chicago,  Dr.  Roger  Gendrin  of  the  French  National  Laboratory  of 
Ionospheric  and  Magnetospheric  Physics,  Professor  James  van  Allen  of  the  University  of 
Iowa,  Professor  Charles  Kennel  of  UCLA  and  the  Scripps  Institution  (and  a  former 
Associate  Administrator  of  NASA),  and  Professor  Alfven  himself  who  inaugurated  the 
series. 

Jointly  with  AFRL,  we  also  organized  a  series  of  workshops  in  "Geoplasma 
Research".  Each  workshop  was  targeted  at  a  specific  topic  of  frontier  geoplasma  physics 
and  includes  basic  tutorial  talks  and  invited  specialty  lectures.  The  format  of  the 
workshops  has  been  designed  to  allow  ample  discussion  time  and  individual  interactions. 
Since  the  untimely  passing  of  Professor  Alfven,  we  have  combined  the  symposium- 
workshop  activities  into  a  single  series  entitled:  the  "Cambridge  Symposia/Workshops  on 
Space  Plasma  Physics".  During  our  recent  15th  Cambridge  Symposium/Workshop  held 
in  Cascais,  Portugal,  Professor  Carl-Gunne  Falthammar  of  the  Royal  Institute  of 
Technology  and  the  Alfven  Laboratory  of  Plasma  Physics,  a  favorite  student  of  Professor 
Alfven,  was  the  first  Alfven  Memorial  Lecturer  of  the  new  series.  Such  activities  have 
received  much  praise  from  the  worldwide  geo/cosmo  community.  Proceedings  of  the 
series  entitled,  "Physics  of  Space  Plasmas",  have  become  informal  textbooks  treasured  by 
both  established  scientists  and  students  in  geo/cosmo  research.  .  The  center  also 
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periodically  co-sponsors  other  scientific  workshops  and  symposia  with  other  institutions 
on  various  important  topics  of  current  space  plasma  research.  One  such  recent  volume 
was  devoted  to  the  memory  of  Professor  S.  Chandrasekhar  of  the  University  of  Chicago. 
Professor  Chandrasekhar,  a  Nobel  laureate,  was  a  giant  in  the  scientific  fields  of  plasma 
physics,  statistical  mechanics,  astrophysics  and  space  sciences,  radiative  transfer,  and 
MHD  stability.  This  memorial  volume  covers  the  entire  spectrum  of  "Chandra’s" 
research  interests.  Other  noteworthy  contributions  by  the  center  include  the  editing  of  the 
collection  of  a  set  of  outstanding  contributions  by  a  group  of  international  scientists 
whose  specialty  are  in  the  physics  of  complexity  in  a  special  issue  entitled  "Forced  and/or 
Self-Organiced  Criticality  [FSOC]  in.  Space  Plasmas"  published  by  the  Journal  of 
Atmospheric  and  Solar  Terrestrial  Physics  [JASTP],  and  a  Special  Topics  Section  on 
Transport  Processes  in  Geospace  in  the  journal  of  Plasma  Physics. 

Our  center  has  interacted  actively  with  a  number  or  research  organizations  and 
universities  in  addition  to  the  Air  Force  Research  Laboratory.  These  include,  the 
International  Space  Science  Institute  in  Bern,  Switzerland,  the  Italian  National  Space 
Research  Institute  in  Rome,  the  Institute  of  Space  and  Astronautical  Science,  Japan,  the 
University  of  Calgary,  Canada,  the  Applied  Physics  Laboratory  of  the  Johns  Hopkins 
University,  the  Max-Planck  Institutes  for  Extraterrestrial  Physics  and  Aeronomy  in 
Germany,  the  Lockheed-Martin  Palo  Alto  Research  Laboratory,  the  Universities  of 
California  at  Berkeley,  Los  Angeles  and  Irvine,  the  Cornell  University,  the  University  of 
New  Hampshire,  NASA  Goddard  Space  Flight  Center,  the  National  Research  Council  of 
Canada,  the  Naval  Research  Laboratory,  the  Imperial  College  of  London,  the  University 
of  Warwick,  England,  the  British  Antarctic  Survey,  Boston  College  and  the  Boston 
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University.  Visits  to  these  institutions  and  by  scientists  from  these  institutions  have 
provided  the  necessary  stimulus  to  keep  our  research  program  vibrant,  up-to-date,  and  at 
the  same  time  constantly  in  touch  with  practical  applications. 

One  of  the  prime  missions  of  the  MIT  Center  for  Theoretical  Geo/Cosmo  Plasma 
Physics  is  to  provide  an  environment  for  the  development  of  young  prospective  students 
in  space  plasma  education.  During  the  past  several  years,  in  cooperation  with  the  MIT 
Graduate  School  and  the  MIT  Minority  Summer  Research  Program  (MSRP),  our  center 
hosted  a  number  of  talented  young  minority  and  female  undergraduate  interns.  By 
intermingling  with  the  established  scientists  at  the  center,  these  young  interns  obtained 
first  hand  knowledge  of  the  true  meaning  of  cutting-edge  scientific  research  in  space 
physics.  The  MIT  Center  also  played  an  active  role  in  the  Research  Scientist  Institute 
(RSI)  program  in  collaboration  with  the  Washington-based  Center  of  Excellence  in 
Education  (CEE)  in  sponsoring  the  education  of  genius-level  high  school  students 
recruited  worldwide  in  preparation  for  them  to  achieve  the  ultimate  of  their  scientific 
capabilities.  CEE  was  founded  by  Admiral  Rickover  in  1983  and  its  honorary  trustees 
include  such  luminaries  as  former  president  Jimmy  Carter,  Senator  Joseph  Leiberman, 
and  General  Colin  Powell.  RSI  alumni  have  dominated  all  scholarship  competitions  of 
national  talent  search  programs  and  many  have  become  leaders  in  the  various  fields  of  the 
21st  century  cutting-edge  scientific  disciplines. 

This  report  is  organized  as  follows:  In  Sec.  Ill,  we  discuss  the  progress  of  our 
center  by  including  sub-sections  of  detailed  descriptions  of  two  of  our  recent  major 
accomplishments,  followed  by  a  complete  list  of  the  center  research  topics.  In  Section 
IV,  we  provide  a  complete  list  of  scientific  publications  by  members  affiliated  with  the 
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Center.  A  complete  list  of  invited  and  review  lectures  by  the  center  personnel  given  at 
various  conferences,  meetings,  and  research  institutions  is  given  in  Section  V.  Finally,  a 
list  of  the  past  and  present  center  personnel  is  given  in  Sec.  VI. 
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III.  PROGRESS  AND  ACCOMPLISHMENTS 

In  our  original  AFOSR  URI  proposal,  we  proposed  a  unique  program  of 
theoretical  research  in  geoplasma  physics.  The  Center  would  be  a  single  cohesive  unit  of 
scientists  from  several  disciplines  interacting  effectively  with  one  another  and  with  groups 
from  external  ongoing  experimental  research  programs.  It  would  not  be  the  purpose  of 
the  Center  to  carry  out  routine  data  analyses.  Instead,  our  approach  would  be  to  interact 
with  the  experimental  groups  and  to  identify  from  the  analyzed  data  those  problems  that 
had  no  ready-made  explanations  and  to  focus  our  efforts  on  the  solution  of  such  new 
problems.  At  all  times,  we  would  not  lose  sight  of  the  practical  applications  of  the 
developed  theories  to  the  prime  mission  of  the  Air  Force. 

During  the  past  few  years  under  the  current  grant  program  “Frontier  Geoplasma 
Research”,  we  have  endeavored  to  follow  such  guidelines  while  developing  the  various 
research  efforts  at  the  Center.  We  identified  many  new,  interesting,  and  at  the  same  time 
puzzling  geoplasma  problems  that  were  not  admissible  to  "routine"  solutions.  We  have 
provided  theoretical  understanding  to  a  number  of  such  identified  problems.  In  many  of 
these  instances,  we  were  able  to  provide  quantitative  descriptions  of  the  phenomena  or 
make  useful  theoretical  predictions  for  future  observations  and  applications. 

The  most  exciting  research  findings  of  our  group  during  this  period  included  the 
development  of  a  first  truly  kinetic  theory  of  the  global  evolution  of  the  fast  stream  solar 
wind  including  the  wave-particle  energization  of  the  ions  and  electrons,  an  "Intermittent 
Turbulence  Theory"  of  sporadic  localized  magnetic  reconnections  that  is  germane  to  the 
coupling  process  between  the  ionosphere  and  magnetosphere  in  the  auroral  region  and  in 
the  dynamics  of  the  plasma  sheet  of  the  Earth's  magnetotail  and  the  solar  wind,  the 
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anomalous  energization  of  ionospheric  ions  in  the  acceleration  and  reverse-current  regions 
of  the  auroral  zone,  and  the  development  of  a  novel  theory  of  the  photoelectron-driven 
polar  wind  based  on  the  kinetic  and  anomalous  heat  transport  of  the  field-aligned 
electrons,  and  nonthermal  ions  in  the  dayside  polar  ionosphere.  We  have  included  below 
two  sub-sections  to  describe  the  details  and  seminal  aspects  of  these  research  results. 

Other  important  research  activities  of  the  center  included:  the  study  of  the 
dynamics  of  the  coherent  structures  such  as  double  layers  and  fast  electron  holes  in  the 
ionosphere,  the  dynamics  of  upward  propagating  lightning  strokes,  the  multi-dimensional 
aspects  of  the  black  auroral  curls,  the  development  and  propagation  of  electromagnetic  ion 
cyclotron  waves  throughout  the  terrestrial  ionosphere,  and  various  processes  related  to  the 
coupling  of  the  ionosphere  and  magnetosphere. 

A  number  of  our  research  findings  have  recently  been  confirmed  by  the 
experimental  data  collected  by  high-time  resolution  instruments  on  board  the  FREJA, 
VIKING,  POLAR,  DE,  FAST  and  "AKEBONO"  satellites  as  well  as  the  TOPAZ, 
SCIFER,  AMICIST,  and  other  high-altitude  rockets.  As  it  has  been  the  history  of  all  past 
research  activities  of  the  MIT  Center,  our  recent  work  is  directly  related  to  AF  C^I 
systems.  The  coupling  phenomenon  between  the  Earth’s  ionosphere  and  magnetosphere, 
intermittent  MHD  turbulence,  the  polar  wind  outflow,  the  black  auroral  curls,  and 
ionospheric  lower  hybrid  cavitary  structures  all  have  profound  effects  on  the  high-latitude 
electron  and  ion  density  profiles  as  well  as  signal  and  wave  propagations.  The 
understanding  of  these  phenomena  is  germane  to  space  surveillance,  space  launch  and 
orbit  operations,  as  well  as  space  weather  modeling  and  forecasting.  In  particular,  our 
understanding  of  pitch-angle  scattering  and  the  dynamics  of  the  central  plasma  sheet 
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allowed  us  to  construct  a  quantitative  model  of  electron  precipitation  in  the  diffuse  aurora. 
Because  this  precipitation  helps  control  the  electron  density  profile  in  the  high-latitude 
ionosphere,  it  has  a  strong  impact  on  the  Air  Force  communication  and  surveillance 
systems  that  must  operate  in  the  region.  Similarly,  because  of  the  effect  of  scintillation  on 
these  systems  and  the  close  relationship  between  scintillations  and  the  high-latitude 
ionospheric  turbulence,  our  quantitative  models  of  the  latter  phenomena  and  its 
consequences  can  be  expected  to  have  great  utility  in  the  practical  business  of  ionospheric 
weather  prediction. 

Listed  below  are  research  topics  that  have  been  studied  and  analyzed  by  the 
members  of  the  center  personnel  during  the  grant  period  followed  by  the  two  subsections 
on  the  solar  wind,  self-organized  criticality,  ion  acceleration,  intermittent  turbulence  and 
complexity. 

•  Anisotropic,  intermittent  turbulence  in  the  solar  wind. 

•  Sporadic,  localized  magnetic  reconnection  and  intermittent  turbulence  in  the 
ionosphere-magnetosphere  coupling  region  and  in  the  Earth's  magnetotail. 

•  Ion  heating  by  low  frequency  waves  in  Earth's  ionosphere  and  magnetosphere. 

•  Kinetic  theory  of  the  global  evolution  of  the  fast  stream  solar  wind  including  the 
wave-particle  energization  of  the  ions  and  electrons. 

•  Lower  hybrid  collapse,  caviton  turbulence,  and  charged  particle  energization  in  the 
topside  ionosphere. 

•  Mode  conversion  processes  involving  the  plasma  turbulence  of  oxygen-hydrogen 
plasmas  in  the  magnetosphere  and  ionosphere. 

•  Anisotropic,  kinetic  polar  wind  driven  by  the  dayside-photoelectrons. 
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Nonlinear  vortex  description  of  black  auroral  curls. 


•  Turbulent  relaxation  of  magnetic  fields  in  space  plasmas. 

•  Ion  and  electron  acceleration  along  auroral  field  lines. 

•  Nonlocal  effects  of  finite  beam-driven  instabilities  in  space  plasmas. 

•  Path  integral  approach  to  nonlinear  particle  acceleration  and  diffusion  in 
ionospheric  plasmas. 

•  Convection  of  ion  cyclotron  waves  to  ion  heating  regions. 

•  Broad-band  spectrum  of  auroral  plasma  turbulence. 

•  Renormalization-group  calculation  of  self-organized  criticality  and  low¬ 
dimensional  behavior  of  auroral  substorm  onsets. 

•  Theory  of  nonlinear  electric  fields  in  the  auroral  acceleration  region. 

•  Stochastic  MHD  models  for  space  plasmas. 

if 

•  Multiple-cyclotron  absorption  of  ion  heating  in  the  cusp/cleft  region. 

•  Wave-particle  ion  cyclotron  turbulence  and  evolution  of  the  electron  distribution  in 
inhomogeneous  space  plasmas. 

•  ULF  waves  along  auroral  field  lines  in  the  central  plasma  sheet. 

•  Energy  source  and  generation  mechanism  for  auroral  kilometric  radiation. 

•  Upward  propagating  intense  lightning  strokes  and  associated  infrared  emissions  in 
the  lower  ionosphere. 

•  Particle  acceleration  by  intense  auroral  VLF  waves. 

•  The  electron  beam  instability  and  turbulence  theories  in  space  plasmas. 

•  Two  stream  interaction  on  auroral  field  lines. 
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•  Energetic  photoelectron  and  the  polar  rain. 

•  Monte-Carlo  modeling  of  polar  wind  electron  distributions  with  anomalous  heat 
flux. 

•  Counterstreaming  electrons  generated  by  lower  hybrid  waves  in  the  auroral  region. 

•  Heating  of  thermal  ions  near  the  equatorward  boundary  of  the  mid-altitude  polar 
cusp. 

•  Stabilization  of  the  cyclotron  autoresonance  maser  instability  in  space  plasmas. 

•  Nonlinear  oblique  whistler  modes  in  collisionless  shocks  of  space  plasmas. 

•  Electromagnetic  tornadoes  in  space. 

•  Radiations  from  large  space  structures  in  low  Earth  orbit  with  induced  AC 
currents. 

•  Ion  waves  and  upgoing  ion  beams  observed  by  the  VIKING  satellite. 

•  Simulation  of  ion  conic  formation  in  the  ionosphere  and  magnetosphere. 

•  Alfven  engine  in  space. 

•  Convection  of  ion  cyclotron  waves  to  ion  heating  regions  in  the  auroral  zone. 

•  Wave  observations  and  their  relation  to  nonresonant  and  resonant  particle  heating 
processes. 
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III.l  Complexity,  Forced  and/or  Self-Organized  Criticality,  and  Topological  Phase 
Transitions  in  Geospace  Plasmas  —  A  Challenge  To  High-Speed  Computing  and 
Numerical  Simulation 
IIL1.1  Introduction 

Research  activity  in  geospace  plasma  physics  is  now  arriving  at  an  interesting 
juncture  that  becomes  apparent  when  looking  back  at  what  has  been  accomplished  in  the 
past  and  looking  forward  to  what  will  be  required  in  the  next  decade.  In  this  regard,  it  is 
noted  that  considerable  observational  and  theoretical  attention  has  been  devoted  towards 
the  understanding  of  local,  point  observations  of  geospace  plasma  phenomena  as  they 
have  been  identified  on  US/European/Japanese  research  satellites.  However,  it  is 
becoming  clear  that  important  questions  that  will  be  receiving  attention  in  the  coming 
years  (particularly  with  the  successful  launches  of  the  CLUSTER  II  and  IMAGE 
satellites)  are  addressed  toward  global  and  multiscale  issues:  questions  of  energy  and 
momentum  transport  within  the  magnetosphere,  ionosphere  and  the  Sun-Earth  connection 
region,  and  of  the  nature  of  the  particle  populations,  their  source,  entry,  energization, 
diffusion,  and  ultimate  loss  from  the  systems.  Thus,  observations  are  becoming  multi- 
spacecraft  and/or  multi-point  in  scope,  and  theoretical  models  are  likewise  being  forced 
to  confront  issues  of  nonlocality  and  "complexity". 

In  this  review,  we  discuss  a  number  of  such  geospace  plasma  processes  of 
complexity,  with  special  emphasis  on  the  dynamics  of  the  magnetotail  and  its  substorm 
behavior.  We  demonstrate  that  the  sporadic  and  localized  interactions  of  magnetic 
coherent  structures  arising  from  plasma  resonances  are  the  origin  of  "complexity"  in 
geospace  plasmas.  Such  interactions,  which  generate  the  anomalous  diffusion,  transport 
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and  evolution  of  the  macroscopic  state  variables  of  the  overall  dynamic  system,  may  be 
modeled  by  a  triggered  localized  chaotic  growth  equation  of  a  set  of  relevant  order 
parameters.  The  dynamics  of  such  intermittent  processes  would  generally  pave  the  way 
for  the  global  system  to  evolve  into  a  "complex"  state  of  long-ranged  interactions  of 
fluctuations,  displaying  the  phenomenon  of  forced  and/or  self-organized  criticality 
(FSOC).  The  coarse-grained  dissipation  due  to  the  intermittent  fluctuations  can  also 
induce  "fluctuation-induced  nonlinear  instabilities"  that  can,  in  turn,  reconfigure  the 
topology  of  background  magnetic  fields. 

The  organization  of  this  paper  is  as  follows.  In  Section  III.  1.2,  we  discuss  the 
origin  of  the  stochastically  distributed  coherent  magnetic  structures  in  the  Earth's 
magnetotail.  In  Section  III.  1.3,  we  consider  the  preferential  acceleration  of  the  coherent 
magnetic  structures  in  the  neutral  sheet  region  and  their  relevance  to  the  interpretation  of 
the  observed  fast  bursty  bulk  flows  (BBF).  In  Section  HI.  1.4,  we  consider  the  localized 
merging  processes  of  the  coherent  structures  and  their  connection  to  the  observed 
localized,  sporadic  reconnection  signatures  in  geoplasmas  (e.g.,  in  the  Earth's 
magnetotail).  The  concept  and  implications  of  the  phenomenon  of  FSOC  are  then 
introduced.  These  ideas  can  lead  naturally  to  the  power  law  scaling  of  the  probability 
distributions  and  fractal  spectra  of  the  intermittent  turbulence  associated  with  the 
coherent  structures.  In  Section  HI.  1.5,  we  provide  some  convincing  arguments  for 
intermittency  and  FSOC  in  the  magnetospheric  plasma  sheet.  In  Section  HI.  1.6,  the 
concept  of  invariant  scaling,  the  dynamic  renormalization-group,  and  topological  phase 
transitions  are  discussed.  These  are  then  related  to  the  multifractal  scaling  across  various 
physical  scales  and  the  fluctuation-induced  nonlinear  instabilities  that  may  lead  to  the 
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triggering  of  global  magnetic  field  reconfigurations  such  as  substorms.  In  Section  III.  1.7, 
simple  phenomenological  and  analogous  intermittency  models  are  described.  The 
summary  and  reference  sections  then  follow. 

III.1.2.  Stochastically  Distributed  Coherent  Structures  in  the  Geospace 
Environment 

In  situ  observations  indicate  that  the  dynamical  processes  in  the  geoplasma 
environment  (e.g.,  the  magnetotail  region  [Angelopoulos  et  al.,  1996;  Lui,  1996;  1998; 
Nagai  et  al.,  1998]),  generally  entail  localized  intermittent  processes  and  anomalous 
global  transports.  It  was  suggested  by  Chang  [1998a, b,c;  1999]  that  instead  of 
considering  the  turbulence  as  a  mixture  of  interacting  waves,  such  type  of  patchy 
intermittency  could  be  more  easily  understood  in  terms  of  the  development,  interaction, 
merging,  preferential  acceleration  and  evolution  of  coherent  magnetic  structures. 

Most  field  theoretical  discussions  begin  with  the  concept  of  propagation  of  waves. 
For  example,  in  the  MHD  formulation,  one  can  combine  the  basic  equations  and  express 
them  in  the  following  propagation  forms: 

pdV/dt=B-VB+-,dB/dt  =  B-Wy+-  (1) 

where  the  ellipsis  represent  the  effects  of  the 
anisotropic  pressure  tensor,  the  compressible 
and  dissipative  effects,  and  all  notations  are 
standard.  Equations  (1)  admit  the  well- 
known  Alfven  waves.  For  such  waves  to 
propagate,  the  propagation  vector  k  must 
contain  a  field-aligned  component,  i.e.. 


Figure  1.  Cross-sectional  view  of  flux  tubes 
in  the  plasma  sheet  of  the  Earth's  magnetotail 
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B-V— >/k-B*0.  However,  at  sites  where  the  parallel  component  of  the  propagation 
vector  vanishes  (the  resonance  sites),  the  fluctuations  are  localized.  That  is,  around  these 
resonance  sites  (usually  in  the  form  of  curves),  it  may  be  shown  that  the  fluctuations  are 
held  back  by  the  background  magnetic  field,  forming  coherent  magnetic  structures  in  the 
form  of  flux  tubes  [Chang,  1998a, b,c;  1999].  For  the 
neutral  sheet  region  of  the  magnetotail,  these  structures 
are  essentially  current  filaments  in  the  cross-tail  direction, 

Fig.  1.  The  results  of  a  2D  MHD  simulation  showing 
such  coherent  structures  are  given  in  Fig.  2. 

Generally,  there  are  various  types  of  propagation 

Figure  2.  2D  simulation 

modes  (whistler  modes,  lower  hybrid  waves,  etc.)  in  a 

geoplasma.  Thus,  we  envision  a  corresponding  number  of  different  types  of  plasma 
resonances  and  associated  coherent  magnetic  structures  that  typically  characterize  the 
dynamics  of  a  plasma  medium  under  the  influence  of  a  background  magnetic  field. 
These  coherent  structures  will  wiggle,  migrate,  deform  and  undergo  different  types  of 
motions  (including  preferential  acceleration),  i.e.,  anisotropic  stochastic  randomization, 
under  the  influence  of  the  local  and  global  plasma  and  magnetic  topology.  It  is  this  type 
of  motion  (evolution)  and  interaction  of  the  coherent  structures  that  characterize  the 
dynamics  of  turbulent  plasma,  not  plane  waves. 

III.1.3.  Preferential  Acceleration  of  BBF 

Consider  the  coherent  magnetic  structures  discussed  in  the  previous  section, 
which  in  the  neutral  sheet  region  of  the  plasma  sheet  are  essentially  filaments  of 
concentrated  currents  in  the  cross-tail  direction.  We  expect  that,  in  a  sheared  magnetic 
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field  that  typically  exists  in  the  region  under  consideration,  these  multiscale  fluctuations 
can  supply  the  required  coarse-graining  dissipation  that  can  produce  nonlinear 
instabilities  leading  to  "X-point-like"  structures  of  the  averaged  magnetic  field  lines.  For 
example,  for  a  sheared  magnetic  field  Bx(z),  upon  the  onset  of  such  fluctuation-induced 
nonlinear  instabilities,  the  averaged  magnetic  field  will  generally  acquire  a  z-component. 
Let  us  choose  x  as  the  Earth-magnetotail  direction  (positive  toward  the  Earth)  and  y  in  the 
cross-tail  current  direction.  Then  the  deformed  magnetic  field  geometry  after  the 
development  of  a  mean  field  "X-point"  magnetic  structure  will  generally  have  a  positive 
(negative)  Bz  component  earthward  (tailward)  of 
the  X-point.  Near  the  neutral  sheet  region,  the 
Lorentz  force  will  therefore  preferentially 
accelerate  the  coherent  structures  (which  are 
essentially  current  filaments  with  polarities 
primarily  pointing  in  the  positive  y-direction) 
earthward  if  they  are  situated  earthward  of  the  X- 
point  and  tailward  if  they  are  situated  tailward  of 
the  X-point,  Fig.  3.  These  results  would  therefore 
match  the  general  directions  of  motion  of  the  observed  BBF  in  the  magnetotail  [Nagai  et 
al,  1998]. 

We  have  performed  preliminary  two-dimensional  numerical  simulations  to  verify 
these  conjectures.  The  simulations  are  based  on  a  compressible  MHD  model  that  has 
been  used  by  us  in  our  previous  studies  of  coherent  magnetic  structures  [Wu  and  Chang, 
2000a, b;  2001]. 


Figure  3.  Schematic  of  a  mean  field  X- 
point  magnetic  field  geometry 
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In  the  first  example,  we  started  the  numerical  simulation  initially  consisting  of  randomly 
distributed  magnetic  fluctuations.  Eventually,  after  some  elapsed  time,  multiscale 


coherent  magnetic  structures  are  formed,  (Figure 
4,  left  top  panel;  left  panels  are  magnetic  fields 
and  right  panels  are  flow  vectors.).  At  this  time, 
a  uniform  Bz  is  introduced.  In  the  middle 
panels  of  Fig.  4,  the  effect  of  the  Lorentz  force  is 
clearly  seen  through  the  motions  of  the  coherent 
structures  after  some  additional  elapsed  time.  If 
the  direction  of  Bz  is  reversed,  the  motions  of 
the  coherent  structures  are  also  reversed  as  seen 
in  the  bottom  panels  of  Fig.  4. 

In  the  second  example,  we  have 


Figure  4.  Effect  of  Lorentz  force  on  the 

motion  of  the  coherent  structures  considered  the  motions  of  the  coherent  magnetic 

structures  that  developed  in  a  sheared  magnetic  field  upon  the  initial  introduction  of 
random  magnetic  fluctuations.  These  structures  were  generally  aligned  in  the  x-direction 
near  the  neutral  line,  z  =  0,  after  some  elapsed 
time.  A  positive  Bz  was  then  applied.  It  can  be 
seen  (from  Fig.  5;  top  panel;  magnetic  fields, 
bottom  panel:  flow  vectors)  that,  after  some 
additional  elapsed  time,  the  coherent  structures 
(mostly  oriented  by  currents  in  the  positive  y- 


Figure  5.  Effect  of  Lorentz  force  on 
direction),  are  generally  accelerated  in  the  positive  coherent  structures  in  a  sheared  field. 
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x-direction;  with  one  exception  where  a  pair  of  magnetic  structures  with  oppositely 
directed  currents  effectively  canceled  out  the  net  effect  of  the  Lorentz  force  on  these 
structures.  We  have  plotted  the  jc-component  of  the  flow  velocities  vx  due  to  the 
cumulative  effect  of  the  Lorentz  (and  pressure)  forces  acting  on  the  flow  (and  in 
particular  the  coherent  structures)  after  some  duration  of  time  has  elapsed  (Fig.  6).  We 


Figure  6.  A  3D  perspective  of  the  vx .  The  peak 
velocity  of  0.02  is  nearly  the  Alfven  velocity. 


note  that  there  are  a  number  of  peaks  and 
valleys  in  the  3-dimensional  display,  with 
peak  velocities  nearly  approaching  that  of 
the  Alfven  speed  mimicking  the  fast  BBF 
that  were  observed  in  the  magnetotail..  It 
is  to  be  noted  that  an  individual  BBF  event 
may  be  composed  of  one,  two,  or  several 
coherent  magnetic  structures. 

In  Fig.7,  we  demonstrate  in  a  self- 
consistent  picture  of  what  might  occur 
near  the  neutral  sheet  region  of  the 
magnetotail.  We  injected  randomly 
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Figure  7.  Contour  plot  of  v,  in  an  X-point  mean 
field  magnetic  geometry. 


distributed  magnetic  fluctuations  in  a 
sheared  magnetic  field.  The  magnetic 
field  geometry  then  underwent  a 
fluctuation-induced  nonlinear  instability 
producing  an  X-point  like  mean  field 
magnetic  structure.  The  coherent 
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structures  (aligned  in  the  neutral  sheet  region)  are  subsequently  accelerated  by  the 
induced  Bz  away  from  the  X-point  in  both  the  positive  and  negative  ^-directions. 
Contour  plots  of  vx  after  some  elapsed  time  clearly  indicate  such  effects. 

III.1.4.  Localized  Reconnections,  FSOC  and  Power  Law  Scaling 

As  the  coherent  magnetic  structures  approach  each  other,  randomly  or  due  to 
external  forcing,  they  may  merge  or  scatter.  Near  the  neutral  sheet  region,  most  of  these 
structures  carry  currents  of  the  same  polarity  (in  the  cross-tail  direction).  As  two  current 
filaments  of  the  same  polarity  migrate  toward  each  other,  a  strong  current  sheet  is 
generated.  The  instabilities  and  turbulence  generated  by  this  strong  current  sheet  can 
then  initiate  the  merging  of  these  structures.  The  results  of  such  merging  processes  might 
be  the  origins  of  the  signatures  of  localized  reconnection  processes  detected  by  ISEE, 
AMPTE,  and  other  spacecraft.  The  observed  localized  reconnection  signatures  to  date 
seem  to  take  place  mainly  in  domain  sizes  comparable  to  that  of  the  ion  gyroradius. 
Thus,  very  probably,  most  of  these  processes  are  influenced  by  microscopic  kinetic 
effects.  During  these  dynamic  processes,  the  ions  can  approximately  be  assumed  to  be 
unmagnetized  and  the  electrons  strongly  magnetized,  and  the  plasma  nearly  collisionless. 
This  can  lead  to  electron-induced  Hall  currents.  With  a  general  magnetic  geometry, 
whistler  fluctuations  may  usually  be  generated. 

Now,  in  analogy  to  the  Alfvenic  resonances,  singularities  of  k\\=  k»B=0  can 
develop  at  which  whistler  fluctuations  cannot  propagate.  These  “whistler  resonances” 
can  then  provide  the  nuclear  sites  for  the  emergence  of  coherent  whistler  magnetic 
structures,  which  are  the  analog  of  the  coherent  Alfvenic  magnetic  structures  but  with 
much  smaller  scale  sizes.  The  intermittent  turbulence  resulting  from  the  intermixing  and 
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interactions  of  the  coherent  whistler  structures  in  the  intense  current  sheet  region  can  then 
provide  the  coarse-grain  averaged  dissipation  that  allows  the  filamentary  current 
structures  to  merge,  interact,  or  breakup.  In  addition  to  the  above  scenarios,  other  plasma 
instabilities  may  set  in  when  conditions  are  favorable  to  initiate  the  merging  and 
interactions. 

FSOC.  As  the  coherent  magnetic  structures  merge  and  evolve,  larger  coherent 
structures  are  formed.  At  the  same  time,  new  fluctuations  of  various  sizes  are  generated. 
These  new  fluctuations  can  provide  the  new  nuclear  sites  for  the  emergence  of  new 
coherent  magnetic  structures.  After  some  elapsed  time,  we  expect  the  distributions  of  the 
sizes  of  the  coherent  structures  to  encompass  nearly  all  observable  scales.  It  has  been 
argued  by  Chang  [1992;  1998a,b,c;  1999]  and  briefly  reviewed  in  Section  6  that,  when 
conditions  are  favorable,  a  state  of  dynamic  criticality  (FSOC)  might  be  approached.  At 
such  a  state,  the  structures  take  on  all  scale  sizes  with  a  power-law  probability 
distribution  of  the  scale  sizes  of  the  fluctuations,  as  well  as  power-law  frequency  (to)  and 
mode  number  (k)  spectra  of  the  correlations  of  the  associated  fluctuations.  Analyses  of 
existing  observations  in  the  intermittent  turbulence  region  of  the  magnetotail  and  those 
conjectured  from  the  AE  index  seem  to  confirm  such  predictions  [Consolini,  1997;  Lui, 
1998;  Angelopoulos  et  al.,  1999;  Lui  et  ah,  2000;  Uritsky  et  ah,  2002b].  In  addition, 
these  multiscale  coherent  structures  may  render  the  coarse-grained  dissipation  that 
sometimes  provides  the  seeds  for  the  excitation  of  nonlinear  instabilities  leading  to,  for 
example,  the  onset  of  substorms  in  the  Earth’s  magnetosphere. 
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III.1.5.  Observational  Evidences  for  FSOC  and  Intermittency  in  Geospace  Plasmas 
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Figure  8.  AE  index  behavior  and  the  LIM  measure. 

The  first  definitive  observation  that  provided  convincing  evidence  indicating 
certain  turbulent  space  plasma  processes  are  in  states  of  forced  and/or  self-organized 
criticality  (FSOC)  was  the  discovery  of  the  apparent  power-law  probability  distribution 
of  solar  flare  intensities  [Lu,  1995].  Recent  statistical  studies  of  complexity  in  geospace 
plasmas  came  from  the  AE  index,  UVI  auroral  imagery,  and  in-situ  measurements  related 
to  the  dynamics  of  the  plasma  sheet  in  the  Earth’s  magnetotail.  For  example,  the  power- 
law  spectra  of  AE  burst  occurrences  as  a  function  of  AE  burst  strength  has  provided  an 
important  indication  that  the  magnetosphere  system  is  generally  in  a  state  of  FSOC 
[Consolini,  1997].  More  recent  evaluations  of  the  local  intermittency  measure  (LIM) 
using  the  Morlett  wavelet  (Figure  8,  from  Consolini  and  Chang  [2001])  clearly  identified 
the  burst  contributions  to  the  AE  index.  These  results  strongly  suggest  that  the  Earth's 
plasma  sheet,  particularly  during  substorms  times,  are  bursty  and  intermittent.  Some  of 
the  salient  features  of  the  complexity  studies  using  the  UVI  imagery  and  in-situ  flow 
measurements  in  the  plasma  sheet  region  verifying  these  suggestions  are  briefly 
described  below. 


Time 
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(i)  Auroral  UVI  Images 

UVI  images  provide  detailed  information  on  the  dynamics  of  spatially  distributed 
magnetotail  activity  covering  extended  observation  periods.  It  has  been  shown  that  the 
positions  of  auroral  active  regions  in  the  nighttime  magnetosphere  are  correlated  with  the 
position  of  the  plasma  sheet  instabilities  [Fairfield  et  al.,  1999;  Lyons  et  al.,  1999; 
Sergeev  et  al.,  1999;  Ieda  et  al.,  2001;  Nakamura  et  al.,  2001a;  2001b],  whereas  timing  of 
the  auroral  disturbances  provides  good  estimates  for  both  small-scale  isolated  plasmoid 
releases  [Ieda  et  al.,  2001]  and  for  the  global-scale  substorm  onset  times  [Germany  et  al., 
1998;  Newell  etal.,  2001]. 

Recently,  Lui  et  al.  [2000]  used  global  auroral  UVI  imagery  from  the  POLAR 
satellite  to  obtain  the  statistics  of  size  and  energy  dissipated  by  .  the  magnetospheric 
system  as  represented  by  the  intensity  of  auroral  emission  on  a  worldwide  scale.  They 
found  that  the  internal  relaxations  of  the  magnetosphere  statistically  follow  power  laws 
that  have  the  same  index  independent  of  the  overall  level  of  activity.  The  analysis 
revealed  two  types  of  energy  dissipation:  those  internal  to  the  magnetosphere  occurring  at 
all  activity  levels  with  no  intrinsic  scale,  and  those  associated  with  active  times 
corresponding  to  global  energy  dissipation  with  a  characteristic  scale. 

Uritsky  et  al.  [2002b]  have  performed  an  extensive  analysis  of  the  probability 
distributions  of  spatiotemporal  magnetospheric  disturbances  as  seen  in  POLAR  UVI 
images  of  the  nighttime  ionosphere. 

They  observed  power-law  statistics  for  the  probability  distributions  of  the 
integrated  sizes  and  energies,  that  are  consistent  with  the  behavior  of  numerical  models 
operating  near  the  FSOC  point;  they  provide  important  observational  evidence  for 


27 


stationary  critical  dynamics  in  the  magnetosphere.  Based  on  this  result,  one  can  expect 
cross-scale  coupling  effects  to  play  a  significant,  if  not  crucial,  role  in  the  development  of 
geomagnetic  disturbances.  More  specifically,  the  FSOC  principle  implies  that  large-scale 
properties  of  the  magnetotail  plasma  sheet  depend  critically  on  the  statistical  hierarchy  of 
small-  and  intermediate-scale  perturbations  associated  with  sporadic  magnetic 
reconnections,  current  sheet  disruptions,  and  other  localized  plasma  instabilities  [Klimas 
et  al.,  2000a,b], 

(ii)  Evidence  of  Intermittency  and  Power-Law  Behavior  of  Bursty  Bulk  Flow  (BBF) 
Durations  in  Earth's  Plasma  Sheet 

Statistical  analyses  based  on  in-situ  data  collected  by  the  GEOTAIL  satellite  have 
yielded  convincing  proof  of  intermittent  turbulence  with  power  law  correlations 
[Angelopoulos  et  al.,  1999].  It  was  shown  that  the  magnetotail  is  generally  in  a  bi-modal 
state:  nearly  stagnant,  except  when  driven  turbulent  by  transport  efficient  fast  flows. 
Figure  9  displays  the  probability  density  functions  of  the  X-component  of  the  plasma 


Figure  9.  Probability  distributions  of  the  non-BBF  and  BBF  x-component  of 
the  normalized  flows  along  with  best  fits  of  Gaussian  and  Castaing  functions. 

sheet  flows  for  both  the  bursty  bulk  flow  (BBF)  and  non-BBF  populations.  Both 

distributions  are  clearly  non-Gaussian  and  both  can  be  fitted  nicely  with  the  Castaing  et 


28 


al.  [1990]  distributions  indicating  intermittency  for  both  components  of  flow  in  the 
plasma  sheet. 


Figure  10.  Probability  density  of 
BBF  bursts  in  the  plasma  sheet. 


Figure  10  shows  the  probability  density  function  of 
the  flow  magnitude  durations  above  400  km/s  in  the 
inner  plasma  sheet  [Angelopoulos  et  al.,  1999].  The 
plasma  sheet  was  selected  on  the  basis  of  plasma 
pressure  ( P.  >  0.01  nPa )  but  very  similar  results  are 

obtained  by  confining  the  data  in  the  near-neutral 
sheet  region  using  plasma  beta  ( >  0.5  ).  A  power 


law  is  clearly  indicated.  The  power-law  behavior  with  spectral  index  ~-1.6  remains  when 
the  velocity  magnitude  threshold  is  changed  to  a  lower  or  higher  value,  and  when 
different  spatial  regions  of  the  magnetotail  plasma  sheet  in  the  X-Y  plane  are  considered, 
indicating  that  the  results  are  quite  robust.  The  results  for  linear  binning  in  flow  duration, 
At,  are  identical  to  the  ones  from  logarithmic  binning  presented  in  the  figure.  Thus, 
there  is  strong  indication  that  the  energy  dissipation  in  the  magnetotail  adheres  to  the 
behavior  expected  from  FSOC. 

III.1.6.  Invariant  Scaling  and  Topological  Phase  Transitions. 

In  the  previous  section,  we  mentioned  the  possibility  of  the  existence  of 
"complex"  topological  states  that  can  exhibit  the  characteristic  phenomenon  of  dynamic 
criticality  similar  to  that  of  equilibrium  phase  transitions.  By  "complex"  topological 
states  we  mean  magnetic  topologies  that  are  not  immediately  deducible  from  the 
elemental  (e.g.,  MHD  and/or  Vlasov)  equations  [Consolini  and  Chang,  2001].  Below,  we 
shall  briefly  address  the  salient  features  of  the  analogy  between  topological  and 
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equilibrium  phase  transitions.  A  thorough  discussion  of  these  ideas  may  be  found  in 
Chang  [1992, 1999;  2001,  and  references  contained  therein.] 

For  nonlinear  stochastic  systems  near  criticality,  the  correlations  among  the 
fluctuations  of  the  random  dynamical  fields  are  extremely  long-ranged  and  there  exist 
many  correlation  scales.  The  dynamics  of  such  systems  are  notoriously  difficult  to 
handle  either  analytically  or  numerically.  On  the  other  hand,  since  the  correlations  are 
extremely  long-ranged,  it  is  reasonable  to  expect  that  the  system  will  exhibit  some  sort  of 
invariance  under  scale  transformations.  A  powerful  technique  that  utilizes  this  invariance 
property  is  the  technique  of  the  dynamic  renormalization-group  [Chang  et  al.,  1978; 
1992;  and  references  contained  therein].  As  it  is  described  in  these  references,  based  on 
the  path  integral  formalism,  the  behavior  of  a  nonlinear  stochastic  system  far  from 
equilibrium  may  be  described  in  terms  of  a  "stochastic  Lagrangian  L".  Then,  the 
renormalization-group  (coarse-graining)  transformation  may  be  formally  expressed  as: 

dL/d£  =  RL  (2) 


where  R  is  the  renormalization-group  (coarse-graining)  transformation  operator  and  t  is 
4  fpj  the  coarse-graining  parameter  for  the  continuous 


group  of  transformation.  It  will  be  convenient  to 
consider  the  state  of  the  stochastic  Lagrangian  in 
terms  of  its  parameters  {P„}.  Equation  (2),  then, 


specifies  how  the  Lagrangian,  L,  flows  (changes)  with 

Figure  11.  Renormalization-group 

trajectories  and  fixed  points.  £  in  the  affine  space  spanned  by  { P„ } ,  Fig.  1 1 . 

Generally,  there  exists  a  number  of  fixed  points  (singular  points)  in  the  flow  field,  where 
dL!di  =  0.  At  a  fixed  point  (L*  or  L**  in  Fig.ll),  the  correlation  length  should  not  be 
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changing.  However,  the  renormalization-group  transformation  requires  that  all  length 
scales  must  change  under  the  coarse-graining  procedure.  Therefore,  to  satisfy  both 
requirements,  the  correlation  length  must  be  either  infinite  or  zero.  When  it  is  at  infinity, 
the  system  is  by  definition  at  criticality.  The  alternative  trivial  case  of  zero  correlation 
length  will  not  be  considered  here. 

To  study  the  stochastic  behavior  of  a  nonlinear  dynamical  system  near  a  particular 
criticality  (e.g.,  the  one  characterized  by  the  fixed  point  L*),  we  can  linearize  the 
renormalization-group  operator  R  about  L*.  The  mathematical  consequence  of  this 
approximation  is  that,  close  to  criticality,  certain  linear  combinations  of  the  parameters 
that  characterize  the  stochastic  Lagrangian  L  will  correlate  with  each  other  in  the  form  of 
power  laws.  This  includes,  in  particular,  the  (k,  0) ),  i.e.  mode  number  and  frequency, 
spectra  of  the  correlations  of  the  various  fluctuations  of  the  dynamic  field  variables.  In 
addition,  it  can  be  demonstrated  from  such  a  linearized  analysis  that  generally  only  a 
small  number  of  (relevant)  parameters  are  needed  to  characterize  the  stochastic  state  of 
the  system  near  criticality  [i.e.,  low-dimensional  behavior;  see  Chang  (1992)]. 

Symmetry  Breaking  and  Topological  Phase  Transitions.  As  the  dynamical 
system  evolves  in  time  (autonomously  or  under  external  forcing),  the  state  of  the  system 
(i.e.,  the  values  of  the  set  of  the  parameters  characterizing  the  stochastic  Lagrangian,  L) 
changes  accordingly.  A  number  of  dynamical  scenarios  are  possible.  For  example,  the 
system  may  evolve  from  a  critical  state  A  (characterized  by  L**)  to  another  critical  state 
B  (characterized  by  L *)  as  shown  in  Fig.  11.  In  this  case,  the  system  may  evolve 
continuously  from  one  critical  state  to  another.  On  the  other  hand,  the  evolution  from  the 
critical  state  A  to  critical  state  C  as  shown  in  Fig.  11  would  probably  involve  a  dynamical 
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instability  characterized  by  a  first-order-like  topological  phase  transition  because  the 
dynamical  path  of  evolution  of  the  stochastic  system  would  have  to  cross  over  a  couple  of 
topological  (renormalization-group)  separatrices.  Alternatively,  a  dynamical  system  may 
evolve  from  a  critical  state  A  to  a  state  D  (as  shown  in  Fig.  11)  which  may  not  be  situated 
in  a  regime  dominated  by  any  of  the  fixed  points;  in  such  a  case,  the  final  state  of  the 
system  will  no  longer  exhibit  any  of  the  characteristic  properties  that  are  associated  with 
dynamic  criticality.  Alternatively,  the  dynamical  system  may  deviate  only  moderately 
from  the  domain  of  a  critical  state  characterized  by  a  particular  fixed  point  such  that  the 
system  may  still  display  low-dimensional  scaling  laws,  but  the  scaling  laws  may  now  be 
deduced  from  straightforward  dimensional  arguments.  The  system  is  then  in  a  so-called 
mean-field  state.  (For  general  references  of  symmetry  breaking  and  nonlinear  crossover, 
see  Chang  and  Stanley  [1973];  Chang  et  al.,  [1973a;  1973b];  Nicoll  et  al.,  [1974;  1976]). 
III.1.7.  Modeling  of  Dynamic  Intermittency 
(i)  Triggered  Local  Order-Disorder  Transitions. 

As  noted  in  the  previous  section,  the  "complex"  behavior  associated  with  the 
intermittent  turbulence  in  magnetized  plasmas  may  generally  be  traced  directly  to  the 
sporadic  and  localized  interactions  of  magnetic  coherent  structures.  Coherent  structures 
may  merge  and  form  more  energetically  favorable  configurations.  They  may  also 
become  unstable  (either  linearly  or  nonlinearly)  and  bifurcate  into  smaller  structures. 
The  origins  of  these  interactions  may  arise  from  the  effects  of  the  various  MHD  and 
kinetic  (linear  and  nonlinear)  instabilities  and  the  resultant  finer-scale  turbulences. 

Such  localized  sporadic  interactions  may  be  modeled  by  the  triggered  (fast) 
localized  chaotic  growth  of  a  set  of  relevant  order  parameters,  0,  (i  =  1, 2 ,N): 
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dOildt  =  yri{OP\cx,c2'--cn\ty,T2”-TH),  (3) 

where  "\jfi "  are  functionals  of  (0,P),  P(x,0  =  Pj(j  =  l,2,..M)  are  the  state  variables  (or 
control  parameters),  (x,t)  are  the  spatial  and  temporal  variables,  the  cn  's  are  a  set  of 
triggering  parameters  and  the  r's  are  the  corresponding  relaxation  time  scales 
characterizing  the  localized  intermittencies,  which,  in  general,  are  much  smaller  than 
those  characterizing  the  evolutionary  time  scales  for  the  state  variables. 

For  the  special  case  of  a  single  one-dimensional  order  parameter  O ,  let  us  assume 
that  the  driving  term  may  be  expressed  in  terms  of  a  real-valued  state  function  (the  local 
configurational  free  energy),  F(0,P;c,,c2  ■■■cn\rvr2  such  that 

dO/dt  =  -dF(0)/dO.  '  (4) 

Generally,  the  topology  of  the  value  of 
the  state  function  F  will  contain  valleys 
and  hills  in  the  real  space  spanned  by  the 
order  parameters  O  at  each  given  (x,t) . 
If  F  is  real,  continuous  with  continuous 
derivatives,  then  we  may  visualize  F 
locally  in  terms  of  a  real  polynomial 
function  of  O.  The  topology  of  the 
state  function  F(0)  may  take  on  forms 
that  are  (1)  locally  stable.  Figure  12a,  (2)  locally  nonlinearly  unstable,  Figure  12b,  and  (3) 
locally  linearly  unstable  (bifurcation).  Figure  12c.  States  (2)  and  (3)  may  be  triggered  by 
a  critical  parameter  c  >  0  (with  a  corresponding  relaxation  time  scale  %  with  or  without 
nonlinear  fluctuations.  For  the  triggering  of  a  locally  nonlinear  instability,  the  triggering 


Figure  12.  Topologies  of  state  function  F(O). 
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parameter  would  typically  be  some  sort  of  a  measure  of  the  amplitude  of  the  local 
nonlinear  fluctuations.  For  the  triggering  of  a  local  linear  instability,  the  triggering 
parameter  would  generally  be  some  sort  of  a  measure  of  the  amplitude  of  the  local 
gradient  of  the  state  variable,  P.  (We  note  that  a  similar  but  symmetrical  model  based  on 
the  Landau-Ginsburg  expansion  has  been  considered  by  Gil  and  Somette  [1996]). 

Alternatively,  a  locally  bifurcated  state  such  as  state  (3)  may  be  influenced  by 
another  triggering  parameter  c'>0  such  that  a  more  preferred  local  state  similar  to  that 
of  state  (1)  becomes  energetically  more  favorable  than  that  of  the  bifurcated  state.  See 
Figure  12d. 

One  may  easily  connect  such  local  triggering  behavior  of  order-disorder 
transitions  to  those  of  the  localized,  sporadic  interactions  of  the  magnetic  coherent 
structures  discussed  in  the  previous  section.  For  example,  the  situation  depicted  by 
Figure  12d  can  be  interpreted  as  the  merging  of  two  coherent  structures.  Also,  states  (2) 
and  (3)  can  be  interpreted  as  the  bifurcation  of  one  coherent  structure  into  two  smaller 
coherent  structures  due  to  certain  nonlinear  or  linear  plasma  instability.  Because  most  of 
such  processes  are  probably  due  to  the  result  of  certain  kinetic  effects,  the  interaction 
time-scales  would  generally  be  of  the  order  of  kinetic  reaction-times  and  therefore  much 
shorter  then  the  system  evolution-time  to  be  described  in  the  following  section. 

(il)  Transport  Equations  of  State  Variables 

The  (slow  time)  evolution  equations  of  the  state  variables  will  generally  consist  of 
convective,  forcing,  random  stirring,  and  transport  terms.  Typically,  they  may  take  on 
the  following  generic  form: 

dPj/dt  =  Cj(P,  0;vvv2,-vm)  (5) 
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where  Cj(j  =  1,2,— M)  are  functionals  of  (P.O)  and  the  v's  are  a  set  of  time  scales 

characterizing  the  long  time  system-wise  evolution  under  the  influence  of  the  chaotic  and 
anomalous  growths  of  the  triggered  order  parameters. 

For  a  one-dimensional  single  real-valued  state  variable,  one  might  envisage  the 
following  typical  transport  equation  in  differential  form: 

dP/dt  =  f  +  h(P,0)dP/dx+Dd2P/dx2  +  g  (6) 

where  /  is  the  forcing  term,  h{P,0)  is  a  convective  function,  g  is  a  random  noise, 
and  D  =  D0H(O)  is  the  anomalous  transport  coefficient  with  H(0),  a  positive-definitive 

growing  function  of  the  magnitude  of  the  absolute  value  of  O . 

From  this  example  (as  well  as  the  general  expression,  (5)),  we  note  that  the 
anomalous  transport  and  convection  can  take  on  varying  magnitudes  (because  of  their 
dependence  on  the  order  parameter)  and  it  is  generally  sporadic  and  localized  throughout 
the  system.  Such  behavior  naturally  leads  to  the  evolution  of  fluctuations  into  all  spatial 
and  temporal  scales,  and  would  generally  lead  the  global  system  to  evolve  into  a 
"complex"  state  of  long-range  interactions  exhibiting  the  phenomenon  of  forced  and/or 
self-organized  criticality  (FSOC)  as  has  been  discussed  extensively  in  several  of  our 
previous  publications  [Chang,  1998a, b,c;  1999;  2001;  Wu  and  Chang,  2000a, b;  2001]. 

(ili)  An  Illustrative  Example. 

We  consider  below  a  simple  2-D  phenomenological  model  to  mimic  non-resonant 
MHD  turbulence.  We  introduce  the  flux  function  yr  for  the  non-resonant  transverse 
fluctuations  such  that  (-3 yf / dx, dy/ 1 dy)  =  (SBy,SBx)  and  V»<5B  =0.  The  coherent 

structures  for  such  a  system  are  generally  flux  tubes  normal  to  the  2-D  plane  such  as 
those  simulated  in  Fig.  2.  Instead  of  invoking  the  standard  2-D  MHD  formalism,  here  we 
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simply  consider  if/  as  a  dynamic  order  parameter.  As  the  flux  tubes  merge  and  interact, 
they  may  correlate  over  long  distances,  which,  in  turn,  will  induce  long  relaxation  times 
near  FSOC.  Assuming  homogeneity,  we  model  the  dynamics  of  flux  tubes,  in  the  crudest 
approximation,  in  terms  of  a  classical  Time-Dependent  Landau-Ginsburg  model  as 
follows: 

dy/kldt=-TkdFldy/_k+fk  (7) 

where  y/k  are  the  Fourier  components  of  the  flux  function,  Fk  an  analytic  function  of 
k2 ,  and  fk  a  random  noise  which  includes  all  the  other  effects  that  had  been  neglected  in 
this  crude  model.  We  shall  assume  the  state  function  to  depend  on  the  flux  function  if/ 
and  the  local  energy  measure  £.  For  non-resonant  fluctuations,  we  shall  assume  that 
diffusion  dominates  over  convection.  Thus,  in  addition  to  the  dynamic  equation  (7),  we 
now  include  a  diffusion  equation  for  % .  In  Fourier  space,  we  have 

Idt  =  -Dk2dFld£_k  +hk  (B) 

where  £k  are  the  Fourier  components  of  £,  D(k)  is  the  diffusion  coefficient, 
F(if/k,£k,k) is  the  state  function,  and  hk  is  a  random  noise.  By  doing  so,  we  separate  the 
slow  energy  transport  due  to  diffusion  of  the  local  energy  measure  £  from  the  noise  term 

fk- 

Under  the  dynamic  renormalization  group  (DRG)  transformation,  the  average  of  the 
square  of  the  fluctuations  of  the  flux  function  in  Fourier  space,  C,  should  scale  as: 

eafc(k,co)  =  C{ket,Q)eaa>e)  (9) 

where  (O  is  the  Fourier  transform  of  the  time  t,  t  the  renormalization  parameter  as 
defined  in  the  previous  section,  and  (ac,aw)  the  correlation  and  dynamic  exponents. 
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Thus,  C/coac/a(0  is  an  absolute  invariant  under  the  DRG,  or  C  ~o)  A ,  where  A  =  -ac  / aw . 
DRG  analyses  performed  for  Gaussian  noises  for  several  approximations  yield  the  value 
for  A  to  be  approximately  equal  to  1.88  -  1.66,  and  a  value  of  -1.0  for  the  (Q  -exponent 
for  the  trace  of  the  transverse  magnetic  correlation  tensor.  Interestingly,  Matthaeus  and 
Goldstein  [1986]  had  suggested  that  such  a  magnetic  correlation  exponent  might 
represent  the  dynamics  of  discrete  structures  in  pseudo  2-D  non-resonant  turbulence; 
thus,  giving  some  credence  to  the  model  and  the  DRG  analysis. 

(Iv)  An  Analogous  Example— the  Global  Modified  Forest-Fire  Model 

While  it  will  be  the  main  goal  of  this  proposed  research  to  develop  physically 
realistic  models  to  study  the  "complexity"  of  the  magnetotail,  we  demonstrate  in  this 
section  an  analogous  model,  which  seems  to  encompass  some  of  the  basic  characteristics 
of  what  are  expected  for  the  dynamical  magnetotail. 

We  describe  below  the  global  generalization  [Tam  et  al.,  2000b]  of  a  "modified 
forest-fire  model"  originally  introduced  by  Bak  et  al.  [1990]  and  modified  by  Drossel  and 
Schwabl  [1992].  Let  us  consider  a  rectangular  grid  of  land,  on  which  trees  may  grow  at 
any  given  site  (i,j)  and  any  time  step  n  with  a  probability  p  .  At  sites  where  there  is  a 
tree,  there  is  a  finite  probability  /  that  it  might  be  hit  by  lightning  and  catch  fire.  If  a 
tree  catches  fire  at  certain  time  step  n ,  then  its  neighbor  will  catch  fire  at  a  subsequent 
time  step  n+ 1 .  It  is  obvious  that  at  any  given  time  step  n ,  there  will  be  patches  of  green 
trees,  patches  of  burning  trees  and  patches  of  empty  spaces.  Fig.  13  (courtesy  of  G. 
Consolini).  One  may  associate  the  growth  of  trees  as  the  development  of  coherent 
structures,  the  burning  trees  as  localized  merging  sites  and  empty  spaces  as  quiescent 
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states  in  the  magnetotail,  a  picture  simulating  the 
sporadically  growing  of  coherent  magnetic  structures 
with  localized  merging. 

The  incremental  changes  of  the  probabilities 
( p,f ),  and  the  densities  of  the  green  trees,  burning 

trees  and  empty  spaces  (A, A, A),  characterize  the 

Figure  13.  Snapshot  of  forest-fire 
simulation. 

dynamics  of  the  forest  fire.  In  particular,  (p,f)  are  the 
dynamic  parameters,  and  the  densities  are  the  dynamic  state  variables.  Following  the 
discussions  given  in  the  previous  section,  we  perform  a  renormalization  (coarse-graining) 
transformation  of  the  parameters  (p,f)  such  that  the  phenomenon  retains  its  essential 
basic  characteristics.  Let  us  symbolically  denote  this  transformation  as: 

p'=F(p,f;pvp2,p3);  f'=G(f,p;pvp2,p3 )  (10) 


When  the  system  reaches  a  steady  state,  the  densities, 
(pvp2,p3) ,  may  be  obtained  from  the  mean  field 

theory  at  any  level  of  coarse-graining  and  are 
expressible  in  terms  of  the  dynamic  parameters  (p,f) . 
In  Figure  14,  we  show  the  affine  space  of 


Figure  14.  Renormalization 
group  trajectories. 


transformations  of  (p,f)  for  a  particular  choice  of 
renormalization  group  procedure  [Tam  et  al.,  2000b]  (a 
global  generalization  of  the  coarse-graining  procedure  suggested  by  Loreto  et  al.  [1995]) 
in  the  physically  meaningful  region  of  (p,f)  between  (0,  1).  We  note  that  there  are  4 
distinct  fixed  points.  Each  fixed  point  has  its  own  distinct  dynamic  critical  behavior.  In 
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the  figure,  the  renormalization  trajectories  are  displayed  as  solid  curves  with  arrows 
indicating  the  direction  of  coarse-graining.  If  the  dynamic  state  of  the  forest-fire  is  at  a 
state  A  near  the  fixed  point  (ii) ,  then  the  system  will  exhibit  the  characteristic  critical 


LogP« 


behavior  prescribed  by  the  invariance  properties  of  this  fixed  point  until  sufficient  coarse- 
graining  has  taken  place  such  that  the  renormalized  system  approaches  a  point  such  as  B 
in  a  region  dominated  by  the  fixed  point  (i) .  At  this  point,  in  the  coarse-grained  view, 
the  dynamical  system  will  now  exhibit  the  critical  behavior  characterized  by  the  fixed 

point  (0 .  We  have  calculated  this  crossover  (symmetry 
breaking)  behavior  for  the  power-law  index  of  the 
probability  distribution  P(s)  of  the  scale  sizes  of  the 
burnt  trees.  Figure  15.  In  terms  of  the  magnetotail 
analogy,  such  type  of  crossover  behavior  might  be 
associated  with  the  change  of  physical  behavior  from 
the  kinetic  state  to  the  MHD  state.  We  note  that  there  is 


Figure  15.  Crossover  between 
two  fixed  points. 


a  distinct  separatrix  connecting  the  fixed  points  (iii)  and  (iv) .  Thus,  transition  from  a 


dynamic  state  such  as  B  to  a  dynamic  state  such  as  C  generally  cannot  be  smooth  and  is 
probably  catastrophic.  This  would  be  the  analog  of  the  triggering  of  the  onset  of 
substorms  due  to  coarse-grained  dissipation  in  the  magnetotail.  Although  not  detailed 
here,  there  are  many  other  properties  of  this  analogous  model  that  resemble  the 
magneotail  dynamics. 
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III.  1.8.  Summary 

We  have  presented  a  dynamical  theory  of  “complexity”  for  space  plasmas  far 
from  equilibrium.  The  theory  is  based  on  the  physical  concepts  of  mutual  interactions  of 
coherent  magnetic  structures  that  emerge  naturally  from  plasma  resonances. 

Models  are  constructed  to  represent  the  local  intermittencies  and  global 
evolutional  processes.  Dynamical  renormalization-group  techniques  are  introduced  to 
handle  the  stochastic  nature  of  the  model  equations  that  exhibit  “complexity”.  Examples 
are  provided  to  illustrate  the  phenomena  of  anomalous  transports,  preferential 
acceleration,  forced  and/or  self-organized  criticality  [FSOC],  symmetry  breaking  and 
general  topological  phase  transitions  including  the  reconfigurations  of  mean  magnetic 
field  geometries  due  to  coarse-grained  dissipation. 

Both  the  physical  concepts  and  mathematical  techniques  are  non  traditional.  The 
readers  are  encouraged  to  consult  the  original  references  for  further  in-depth  studies  of 
these  new  emerging  ideas. 
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III.2  Comparison  of  the  effects  of  wave-particle  interactions  and  the  ki¬ 
netic  suprathermal  electron  population  on  the  acceleration  of  the 
solar  wind 

III.2.1  Introduction 

Since  the  pioneering  work  by  Belcher  and  Davis  [1971],  the  existence  of  waves  has 
been  well  established  in  the  research  of  the  fast  solar  wind.  Wave-particle  interactions 
have  been  considered  as  a  possible  solar  wind  heating  mechanism,  and  their  effects 
have  been  studied  by  a  number  of  authors  [e.g.  Hollweg  and  Turner,  1978;  Dusenbery 
and  Hollweg,  1981;  Marsch  et  al.,  1982c;  Isenberg,  1984;  Leer  et  al.,  1982;  Hu  et  al., 
1997;  Czechowski  et  al.,  1998;  Cranmer  et  al.,  1999b;  Hu  and  Habbal,  1999;  Li  et 
al.,  1999,  Tam  and  Chang,  1999a;  Isenberg  et  al.,  2000;  Cranmer,  2001;  Gary  et  al., 
2001;  Liewer  et  al.,  2001;  Tam  and  Chang,  2001;  Tu  and  Marsch,  2001;  Vocks  and 
Marsch,  2001].  Experimentally,  there  has  been  increasing  evidence  suggesting  that 
wave-particle  interactions  play  a  significant  role  in  the  evolution  and  acceleration 
of  the  solar  wind.  First,  the  proton  velocity  distributions  observed  by  the  Helios 
spacecraft  featured  an  anisotropy  in  the  core  region,  where  the  temperature  in  the 
perpendicular  direction  was  higher  than  that  in  the  direction  parallel  to  the  magnetic 
field  [ Marsch  et  al,  1982b].  Such  an  anisotropy  was  an  indication  of  the  presence  of  ion 
perpendicular  heating  in  the  solar  wind.  The  observed  proton  velocity  distributions 
were  also  frequently  characterized  by  double  peaks  in  the  field-aligned  direction.  The 
formation  of  these  double  peaks  has  recently  been  associated  with  kinetic  ion  cyclotron 
resonances  involving  sunward  propagating  electromagnetic  waves  [Tam  and  Chang , 
1999a, b]. 

The  differential  speeds  between  the  alpha  particles  and  the  protons  in  the  solar 
wind,  which  were  perhaps  the  most  unexpected  experimental  results  observed  by  the 
Helios  spacecraft,  may  also  be  highly  related  to  wave-particle  interactions.  Measure- 
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ments  by  the  spacecraft  revealed  that  the  helium  ions  have  a  higher  outflow  speed 
compared  with  the  protons  [Marsch  et  al. ,  1982a],  especially  in  the  fast  solar  wind. 
Indeed,  more  recent  observations  by  the  WIND  and  the  Ulysses  spacecraft  have  also 
arrived  at  a  similar  conclusion  [Steinberg  et  al.,  1996;  Feldman  et  al.,  1996].  Because 
the  helium  ions  have  a  smaller  charge-to-mass  ratio  compared  with  the  protons,  the 
observed  ion  differential  speeds  were  inconsistent  with  the  idea  of  the  field-aligned 
ambipolar  electric  field  being  the  main  solar  wind  acceleration  mechanism.  A  phys¬ 
ical  mechanism  that  has  been  demonstrated  to  preferentially  accelerate  helium  ions 
rather  than  protons  is  kinetic  ion  resonant  heating  by  electromagnetic  waves  in  the 
ion  gyrofrequency  range  [Tam  and  Chang ,  1999a, b].  These  waves  may  be  produced 
by  the  small-scale  reconnections  near  the  coronal  region.  Recently,  Chang  [2003]  has 
suggested  that  such  type  of  intermittent  turbulence  may  be  associated  with  the  “com¬ 
plexity”  generated  by  the  sporadic  localized  mergings  and  interactions  of  the  coherent 
magnetic  structures  near  the  coronal  holes.  Measurements  of  t.he  power  spectra  of 
these  waves  by  the  Helios  spacecraft  indicated  that  the  wave  power  decreases  with 
frequency  [Bavassano  et  al.,  1982;  Denskat  and  Neubauer,  1982].  Note  that  the  alpha 
particles  have  a  lower  gyrofrequency  than  the  protons.  Therefore,  in  general,  more 
wave  power  is  available  to  the  helium  ions  than  to  the  protons,  at  least  in  the  ob¬ 
served  regimes  (0.3-1  AU).  Extrapolating  the  power  spectrum  to  distances  close  to 
the  sun,  Tam  and  Chang  [1999a]  have  shown  that  the  kinetic  wave-particle  interac¬ 
tions  in  general  transfer  more  energy  to  a  helium  ion  than  to  a  proton.  Although 
the  energy  is  mainly  transferred  into  the  transverse  direction,  it  converts  into  the 
field-aligned  direction  downstream  in  the  solar  wind  due  to  the  mirror  effect,  leading 
to  a  preferential  acceleration  of  the  helium  ions. 

Experimental  data  from  SOHO  have  provided  even  more  evidence  in  support  of 
the  significance  of  wave-particle  interactions  in  the  fast  solar  wind.  The  observed  05+ 
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ions  in  the  extended  corona  deviate  significantly  from  those  of  a  thermal  distribution. 
Their  velocity  distributions  are  strongly  anisotropic,  with  T±  >  Tj|,  and  perpendicular 
kinetic  temperature  higher  than  108  K  [Kohl  et  al.,  1998;  Cranmer  et  al .,  1998,  1999a]. 
It  has  been  pointed  out  by  the  investigators  who  analyzed  the  data  that  this  evidence 
is  consistent  with  ion  cyclotron  resonant  heating. 

Besides  wave-particle  interactions,  other  physical  mechanisms  have  also  been  con¬ 
sidered  as  possible  drivers  of  the  solar  wind.  The  proposal  of  suprathermal  electron 
effects  has  led  to  theory-data  agreement  in  the  ionospheric  polar  wind,  which,  like 
the  solar  wind,  is  an  outflow  of  space  plasmas  along  open  magnetic  field  lines.  The 
reason  why  suprathermal  electrons  may  drive  an  outflowing  plasma,  such  as  the  solar 
wind,  is  the  following.  Without  taking  wave-particle  interactions  into  consideration, 
Scudder  and  Olbert  [1979]  argued  that  due  to  the  velocity-dependence  of  the  Coulomb 
collisional  depth  and  the  global  kinetic  nature  of  the  solar  wind  outflow,  suprather¬ 
mal  tails  could  be  formed  in  the  electron  distributions.  To  demonstrate  this  “velocity 
filtration  effect”,  these  authors  used  a  simplified  collisional  operator  to  determine  the 
evolution  of  an  electron  distribution  along  the  solar  wind.  They  were  able  to  show 
that  an  electron  distribution  in  the  form  of  a  Maxwellian  at  the  corona  could  evolve 
into  one  with  a  notable  tail  population  in  the  anti-sunward  direction  near  1  AU.  Olbert 
[1982]  further  suggested  that  the  heat  flux  contribution  by  the  suprathermal  electron 
tail  may  drive  the  solar  wind.  Such  an  idea  has  been  applied  to  the  ionospheric  po¬ 
lar  wind  with  photoelectrons  playing  the  role  of  the  suprathermal  population  [Tam 
et  al.,  1995,  1998],  and  has  successfully  addressed  the  various  features  observed  by 
satellites.  It  has  been  shown  that  in  the  absence  of  wave-particle  interactions,  the 
suprathermal  electron  population  can  increase  the  ambipolar  electric  field,  leading  to 
higher  ion  outflow  speeds  in  the  polar  wind.  Because  the  ionospheric  polar  wind  and 
the  solar  wind  are  analogous  phenomena  —  both  are  outflow  of  space  plasmas  along 
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open  magnetic  field  lines  —  it  is  reasonable  to  expect  that  suprathermal  electron 
effects  may  accelerate  the  ions  in  the  solar  wind. 

The  combined  effects  of  ion  cyclotron  resonant  heating  and  the  suprathermal 
electron  population  have  been  considered  in  a  solar  wind  study  by  Tam  and  Chang 
[1999a].  The  study  was  the  first  to  describe  the  global  kinetic  evolution  and  accel¬ 
eration  of  the  fast  solar  wind  under  the  influence  of  the  wave-particle  interactions, 
Coulomb  collisions,  and  an  ambipolar  electric  field  that  was  consistent  with  the  parti¬ 
cle  distributions  themselves.  Results  of  the  investigation  associated  the  wave-particle 
interactions  with  some  of  the  solar  wind  observations  discussed  earlier,  including  the 
differential  speed  between  the  helium  ions  and  the  protons,  and  the  double-peaked 
proton  velocity  distributions.  The  bulk  acceleration  of  the  solar  wind  demonstrated 
in  the  study,  however,  consisted  of  contributions  due  to  both  the  ion  resonances  and 
the  suprathermal  electron  effects.  It  was  therefore  unclear  to  what  extent  each  of 
the  two  acceleration  mechanisms  contributes  to  the  driving  of  the  solar  wind.  Our 
goal  in  this  study  is  to  compare  the  relative  importance  of  the  kinetic  effects  due  to 
suprathermal  electrons  and  ion  cyclotron  resonances,  and  to  determine  which  of  the 
acceleration  mechanisms  is  mostly  responsible  for  the  driving  of  the  solar  wind. 

III. 2. 2  Model 

Our  model  in  this  study  is  based  on  a  self-consistent  hybrid  technique  originally 
developed  to  describe  the  ionospheric  polar  wind  [ Tam  et  al.,  1995,  1998].  Because 
the  polar  wind  (as  termed  in  analogy  with  the  solar  wind  [ Axford ,  1968;  Banks  and 
Holzer ,  1968])  is  also  an  example  of  plasma  outflow  along  open  magnetic  field  lines, 
our  present  model  retains  the  general  features  of  the  self-consistent  hybrid  model: 
It  is  based  on  an  iterative  scheme  between  fluid  and  kinetic  calculations.  In  the 
fluid  calculations,  a  set  of  equations  determines  the  ambipolar  electric  field  and  the 
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properties  of  the  bulk  thermal  electrons,  whose  distributions  are  assumed  to  be  in 
a  drifting  Maxwellian.  To  solve  the  fluid  equations,  we  impose  the  quasi-neutrality 
and  current-free  constraints  for  the  solar  wind,  and  make  use  of  the  results  from  the 
kinetic  calculations  in  the  model.  The  kinetic  part  of  the  model  consists  of  multiple 
components,  each  describing  the  evolution  of  a  particle  component  along  the  solar 
wind.  In  the  basic  version  of  the  model,  the  kinetic  approach  is  applied  only  to  two 
particle  components:  proton  and  alpha  particle,  the  two  major  ion  species  in  the 
solar  wind.  The  distributions  of  these  ions  evolve  under  the  influence  of  Coulomb 
collisions,  an  ambipolar  electric  field,  and  wave-particle  interactions.  Results  based 
on  the  global  evolution  of  these  distributions  are  coupled  into  the  fluid  calculations. 
This  coupling  ensures  the  consistency  of  the  ambipolar  electric  field  with  the  particle 
distributions  when  an  iteration  between  the  fluid  and  kinetic  parts  of  the  model 
converges.  In  addition,  a  convergence  of  the  results  enables  our  kinetic  calculations 
to  correctly  take  into  account  the  Coulomb  collisions  for  the  ions,  including  those 
among  the  same  ion  species. 

Effects  of  kinetic  ion  resonant  heating  are  incorporated  into  the  model  with  a 
Monte  Carlo  technique.  This  technique  was  originally  developed  to  model  ion  resonant 
heating  in  the  auroral  region,  and  has  successfully  addressed  the  existence  of  ion 
conics  [Retterer  et  al.,  1983].  By  including  this  technique  into  the  model,  we  are 
able  to  incorporate  the  resonant  effect  of  a  given  wave  spectrum  into  an  otherwise 
self-consistent  solar  wind  description. 

To  be  more  quantitative  on  the  description  of  our  kinetic  calculations,  we  solve 
the  following  steady-state  collisional  kinetic  equations  for  the  protons  and  the  helium 
ions,  the  particle  components  that  resonate  with  the  waves: 
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where  s  is  the  distance  along  the  radial  magnetic  field  line,  vj|  and  v±_  denote  velocities 
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in  the  spacecraft  frame,  fj(s,  U||,  Vx)  is  the  distribution  function  for  the  particle  species 
j,  q  and  m  are  the  electric  charge  and  mass  of  the  species  respectively,  E\\  is  the  field- 
aligned  ambipolar  electric  field,  g  is  the  gravitational  acceleration,  B  is  the  magnetic 
field,  B'  =  dB/ds ,  Cj  is  a  Coulomb  collisional  operator  for  the  species  j,  and 
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(III. 2. 2) 


is  an  operator  that  describes  resonant  heating,  with  Dj\\  and  Djx  being  quasilinear 
diffusion  coefficients  that  represent  heating  in  the  parallel  and  perpendicular  direc¬ 
tions  respectively.  The  expressions  for  Dj\\  and  Djx ,  in  general,  depend  on  a  number 
of  factors,  including  the  kind  of  wave-particle  interactions,  the  conditions  for  reso¬ 
nances,  the  power  of  the  waves,  etc.  Below,  we  discuss  in  detail  the  conditions  that 
determine  the  diffusion  coefficients  in  the  study. 

Ion  Cyclotron  Resonant  Heating 

The  wave-particle  interactions  to  be  described  in  Eq.  (III.2.1)  and  (III.2.2)  are 
kinetic  ion  cyclotron  resonances,  which  have  been  demonstrated  to  preferentially 
accelerate  the  helium  ions  over  the  protons  [  Tam  and  Chang ,  1999a, b].  The  reso¬ 
nances  that  we  consider  involve  electromagnetic  waves  whose  wavevector  component 
kL  <C  fc||  ft  k.  Because  the  thermal  speeds  of  the  ions  are  much  smaller  than  the  local 
Alfven  speed  near  the  corona,  where  the  majority  of  the  resonant  heating  occurs,  the 
properties  of  the  waves  can  be  approximated  by  a  cold  plasma  dispersion  relation. 
Moreover,  due  to  the  relatively  low  thermal  speeds,  few  ions  are  able  to  resonate  with 
right-hand  polarized  waves.  Therefore,  the  only  waves  that  we  consider  in  this  study 
are  left-hand  polarized,  and  are  in  the  Alfven/ion  cyclotron  branch,  which  in  fact  con¬ 
sists  of  two  sub-branches  for  the  two-ion  (proton  and  helium)  plasma  in  our  model. 
In  the  solar  wind  frame,  the  lower  sub-branch  covers  all  frequencies  up  to  the  helium 
gyrofrequency,  f la.  The  upper  sub-branch  plateaus  at  the  proton  gyrofrequency  f lp 
for  large  k ,  and  its  frequency  limit  at  infinite  wavelength  is  above  fla,  leaving  a  gap 
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between  the  two  sub-branches.  For  an  ion  (proton  or  alpha  particle)  with  parallel 
velocity  v'\\  in  the  solar  wind  frame  to  resonate  with  the  left-hand  polarized  waves, 
the  following  condition  has  to  be  satisfied: 

J  -  fcv'n  =  0,  (III. 2. 3) 

where  u'  is  the  wave  frequency  in  the  solar  wind  frame,  and  f I  is  the  gyrofrequency 
of  the  ion.  The  wave  frequency  observed  by  the  spacecraft  is  Doppler-shifted  to  u, 
given  by: 

u  =  J  +  kU,  (III. 2.4) 

where  U  is  the  local  solar  wind  speed.  Note  that  the  sign  of  k  can  be  positive 
or  negative,  representing  outward  (anti-sunward)  or  inward  (sunward)  propagating 
waves  respectively. 

Although  the  waves  in  the  solar  wind  propagate  predominantly  outward,  a  sta¬ 
tistical  analysis  of  the  power  spectra  measured  by  the  Helios  spacecraft  indicated 
that  inward  propagating  waves  constitute  a  considerable  amount  of  power  in  both 
the  fast  and  slow  solar  wind  [Tu  and  Marsch,  1990].  The  study  showed  that  in  the 
fast  solar  wind,  the  power  of  inward  propagating  waves  (e_)  tends  to  approach  its 
outward  counterpart  (e+)  as  the  radial  distance  increases.  The  results  of  the  study 
motivated  us  to  include  inward  propagating  waves  in  our  model.  We  assume  the 
following  relation  for  the  wave  power  ratio  in  the  interplanetary  space: 

e_/e+  =  1  -  e-(s/3>2  ,  (III.2.5) 

where  s  is  the  heliocentric  distance  in  AU.  Note  that  with  Eq.  (III. 2.5),  we  assume 
that  the  inward  propagating  waves  constitute  only  2.4  x  10-6  of  the  total  wave  power 
at  1  R©  (solar  radius).  Even  with  such  a  small  ratio  of  e_/e+  near  the  corona,  the 
presence  of  these  sunward  propagating  waves  have  been  shown  to  lead  to  double- 
peaked  proton  velocity  distributions  in  the  solar  wind  [Tam  and  Chang ,  1999a].  The 

reason  for  the  formation  of  the  non-thermal  feature  is  as  follows.  Because  u'  <  Qp  for 
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all  the  waves  in  the  ion  cyclotron  branch,  it  is  clear  from  Eq.  (III. 2. 3)  that  protons 
with  v'\\  <  0  can  only  resonate  with  the  outward  propagating  waves,  while  those  with 
t/||  >  0  only  interact  with  the  incoming  waves.  Note  that  at  any  given  distance,  the 
anti-sunward  and  sunward  propagating  waves  heat  two  distinct  proton  populations 
in  the  velocity  space.  If  incoming  waves  were  absent,  pitch-angle  diffusion  of  the 
protons  due  to  the  wave-particle  interactions  could  only  occur  at  velocity  space  where 
i/||  <  0.  The  presence  of  sunward  propagating  waves  enables  the  proton  distributions 
to  undergo  two  distinct  sets  of  pitch-angle  diffusion,  one  at  v'\\  <  0  and  the  other  at 
v'\\  >  0.  As  the  solar  wind  evolves  along  the  flow,  each  of  these  two  sets  of  velocity- 
space  diffusion  gives  rise  to  a  hump  in  the  proton  distributions  because  of  the  effect  of 
mirror  folding,  thereby  forming  the  observed  double-peaked  feature.  Later,  we  shall 
provide  evidence  to  support  the  association  between  double-peaked  proton  velocity 
distributions  and  sunward  propagating  waves  when  we  briefly  discuss  the  effect  on 
the  solar  wind  due  to  a  different  e_/e+  ratio. 

The  power  of  the  waves  in  the  solar  wind  decreases  with  frequency,  as  indicated 
by  the  Helios  measurements  of  the  magnetic  field  spectra  [Bavassano  et  ai,  1982; 
Denskat  and  Neubauer,  1982].  The  dependence  of  the  wave  power  on  frequency  can 
generally  be  described  by  a  power  law.  It  was  also  found  that  the  wave  power  decreases 
with  radial  distance,  and  the  dependence  can  also  be  described  by  a  power  law,  at 
least  in  the  observed  regime  [ Bavassano  et  al. ,  1982].  In  our  model,  we  assume  a 
magnetic  field  power  spectrum  based  on  these  experimental  results.  We  extend  the 
power-law  relation  to  higher  frequencies.  As  for  the  dependence  on  radial  distance, 
we  assume  a  power-law  relation  and  interpolate  to  obtain  the  magnetic  field  power 
between  the  two  experimental  limits  (0.29  and  0.87  AU).  For  radial  distances  outside 
the  measurement  range,  we  assume  the  wave  power  to  be  proportional  to  s~4'2. 

With  the  resonance  condition  and  magnetic  field  power  spectra  defined,  we  are 
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able  to  determine  the  frequencies  and  power  of  the  waves  that  are  in  resonance  with 
the  ions.  The  two-ion  (proton  and  helium)  cold  plasma  dispersion  relation  for  the 
left-hand  polarized  waves  in  the  solar  wind  frame,  which  we  consider  to  be  the  frame 
where  both  the  combined  ion  populations  and  the  total  electron  population  carry  zero 
current,  can  be  expressed  as: 


(■ kvA f  =  u)'2 


1  —  u)1  juj 0 

(1  -wVO,)(l  --(S/Slp) J  ’ 


where  va  is  the  Alfven  speed,  and 


(III. 2. 6) 


loq  —  Oo.(np  +  4 T  2uq, ) 


(III. 2. 7) 


is  a  frequency  that  depends  on  the  densities  of  the  protons  (np)  and  the  helium 
ions  (na).  Note  that  u>0  represents  the  long- wavelength  limit  for  the  upper  ion  cy¬ 
clotron  sub-branch.  One  can  see  that  Eq.  (III.2.6)  does  not  have  any  solution  for 
£la  <  u>'  <  u>o,  resulting  in  a  frequency  gap  between  the  two  ion  cyclotron  sub¬ 
branches.  However,  in  the  limit  of  a  one-ion  plasma  (na  — >■  0),  the  frequency  gap 
disappears  (as  u>0  -4  Cla  in  Eq.  (III. 2. 7)),  and  the  dispersion  relation  reduces  to  a 
single  branch.  Equation  (III.2.6),  together  with  Eq.  (III.2.3),  determine  the  wavevec- 
tors  and  frequencies  of  the  waves  that  resonate  with  the  individual  ions  in  the  solar 
wind  frame.  These  resonant  frequencies  are  then  Doppler-shifted  to  the  spacecraft 
frame  according  to  Eq.  (III. 2.4).  The  total  wave  power  at  these  frequencies  are  then 
determined  based  on  the  interpolation/extrapolation  scheme  of  the  magnetic  field 
spectrum  described  earlier. 

The  diffusion  coefficients  in  Eq.  (III. 2.2)  involve  all  of  the  considerations  discussed 
above.  Standard  expressions  for  these  coefficients  can  be  found  in,  for  example, 
Swanson  [1989].  But  by  taking  into  account  all  the  considerations  above,  we  arrive 
at  a  modified  version  of  the  diffusion  coefficients: 

■Dii  =  n  ^  [  duvA2PB(uj/2'iT)5(u!  —  kv\i  — Q,),  (III. 2. 8) 

'  Am,2  J 
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Dj_  =  f  du  {n/k)2PB{u/2n)  5(u  -  kv N  -  0),  (III.2.9) 

where  the  (^-functions  reflect  our  emphasis  on  the  resonance  condition  equivalent  to 
Eq.  (III. 2. 3),  u>  and  k  are  related  through  Eq.  (III. 2. 3),  (III. 2.4)  and  (III. 2. 6),  PB 
is  a  function  for  the  magnetic  field  wave  power,  which  is  based  on  our  interpola¬ 
tion/extrapolation  scheme  of  the  Helios  measurements,  and  depends  on  the  sign  of  k 
(direction  of  the  wave  propagation)  as  described  earlier,  and  lastly, 

77  =  770  exp  [(1  —  r)/0.5]  (III. 2. 10) 

is  an  efficiency  factor  that  adjusts  the  available  wave  power,  and  depends  on  the 
heliocentric  distance  r  (in  the  unit  of  R0).  The  adjustment  of  the  wave  power  is 
mainly  for  two  reasons.  First,  the  factor  rjo  in  Eq.  (III. 2. 10)  accounts  for  the  fact  that 
the  left-hand  circularly  polarized  waves  constitute  only  a  fraction  of  the  total  observed 
wave  power.  Second,  the  remaining  factor  on  the  right-hand  side  of  the  equation 
approximates  the  dissipation  of  the  waves  as  they  propagate  along  the  magnetic  field 
line  from  the  sun.  It  is  chosen  based  on  the  study  by  Cranmer  et  al.  [1999b],  who 
suggested  that  most  of  the  waves  originating  from  the  solar  surface  are  dissipated 
before  reaching  1.5  R0. 

Kinetic  Suprathermal  Electron  Effects 

Up  to  this  point,  we  have  described  the  basic  constituents  of  our  model.  This 
model  enables  us  to  follow  the  global  evolution  of  the  ion  distributions  along  the 
solar  wind  flow,  while  taking  into  account  the  influence  due  to  kinetic  ion  cyclotron 
resonant  heating,  Coulomb  collisions,  and  an  ambipolar  electric  field  that  is  consistent 
with  the  particle  distributions  themselves.  The  kinetic  effects  of  the  ions  are  also 
incorporated  into  the  overall  solar  wind  picture  through  the  coupling  of  their  results 
with  the  fluid  calculations.  In  fact,  our  modeling  technique  can  readily  incorporate 
the  kinetic  effects  of  another  particle  component.  To  do  so,  we  simply  add  another 
kinetic  component  into  the  model  to  describe  the  particle  population,  and  modify  the 

60 


quasi-neutrality  and  current-free  conditions,  as  well  as  the  fluid  equations  accordingly. 
This  flexibility  makes  our  model  an  ideal  tool  to  investigate  the  global  kinetic  effects 
of  suprathermal  electrons  on. the  solar  wind. 

To  take  into  account  the  global  kinetic  effects  of  the  suprathermal  electrons,  we 
solve  a  steady-state  collisional  kinetic  equation  for  the  population  that  comprises  the 
tail  portion  of  the  electron  Maxwellian  distribution  at  the  lower  boundary  of  our 
model.  The  equation  is  similar  to  Eq.  (III.2.1),  but  without  the  term  that  represents 
wave-particle  interactions.  In  addition,  the  operator  that  describes  the  Coulomb 
collisions  for  the  suprathermal  electrons  is  simplified.  The  collisional  operator  in 
Eq.  (III.2.1),  such  as  the  Fokker-Planck  collisional  operator,  should  have  the  following 
expression  in  its  complete  form  [see,  e.g.,  Ichimaru,  1973;  Book ,  1989]: 


C,/,(v)  =  -£VvJjb,  (III-2.il) 

b 

where  the  subscript  b  represents  the  background  species  with  which  the  particles  of 


species  j  interact,  and 


Jjb(v)  =  2*Ait  ^ffd3* 
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(III. 2. 12) 

where  A jb  is  the  Coulomb  logarithm,  and  I  is  the  unit  dyad.  In  the  full  collisional 
description,  such  as  that  for  the  ions,  the  background  includes  the  thermal  electrons 
and  the  two  major  ion  species.  But  for  the  suprathermal  electrons,  we  take  into 
account  their  Coulomb  interaction  only  with  the  thermal  electrons.  The  interaction 
among  the  suprathermal  electrons  themselves  are  not  considered  due  to  their  low 
relative  density.  With  these  approximations,  and  the  assumption  that  the  background 
thermal  electron  distribution  is  in  the  form  of  a  drifting  Maxwellian, 

/.(V)  =  o.fe'Tex P  [-m'(V~  “--)-1  ,  (III.2.13) 


the  collisional  term  in  the  kinetic  equation  for  the  suprathermal  electrons  becomes: 


Csfs  —  — Vv  •  (&coii,s/s)  "H  •  (Dcoll,S  •  Vv/s),  (III. 2. 14) 
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where 
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W o 
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aCoii,s  =  -{ms/me)i/0c'f!(xe), 

*M  +  |(c2i_cc)[(i -£)*(..) +  *(4 


with 


(111. 2. 15) 

(111.2.16) 


c  =  V  -  Ue, 

(III.2.17) 

i/0  =  4TrAseneqs2qe2/ms2c3, 

(III.2.18) 

xe  =  mec2 /2Te, 

(III.2.19) 

\P(®)  =  -4=  /  dty/te-*, 
y/n  Jo 

(III.2.20) 

=  dW/dx,  and  the  subscripts  s  and  e  stand  for  the  suprathermal  and  thermal 
electrons  respectively  (thus,  qs  =  qe  and  ms  =  me ).  Equation  (III. 2. 14)  indicates  that 
the  Coulomb  collisions  have  two  effects  on  the  suprathermal  electrons:  to  provide  a  de¬ 
celeration  and  a  diffusion.  The  expression  for  acoii,s  in  Eq.  (III. 2.15)  suggests  that  the 
collisions  act  as  a  friction  between  the  thermal  electrons  and  the  individual  suprather¬ 
mal  electrons.  The  Coulomb  interaction  also  leads  to  a  diffusion  of  the  suprathermal 
electrons  in  the  velocity  space,  characterized  by  the  tensor  Dcoii,s-  Note  that  mathe¬ 
matically,  DcohiS  plays  a  similar  role  as  Dj\\  and  Dj±  of  Eq.  (III.2.2)  in  the  governing 
kinetic  equation.  Therefore,  the  Monte  Carlo  technique  in  our  model  can  also  describe 
the  diffusion  due  to  the  Coulomb  collisions.  According  to  Eq.  (III. 2. 15)  and  (III.2.16), 
both  the  diffusion  and  the  friction  are  characterized  by  a  collisional  frequency  Uq.  Be¬ 
cause  i/0  is  proportional  to  1/c3  (see  Eq.  (III. 2.18)),  the  more  energetic  suprathermal 
electrons  are  less  collisional.  Because  of  this  relationship,  Coulomb  collisions  are  ex¬ 
pected  to  enhance  the  tail  portion  of  the  total  electron  distribution  as  the  particles 
travel  downstream  along  the  solar  wind  flow.  In  our  model,  the  electrons  at  the  tail 
portion  of  the  distribution  have  to  satisfy  two  criteria  at  1  R©,  the  lower  boundary 
in  the  calculations,  in  order  to  be  considered  suprathermal:  (l/2)meu2  >  8Teo  and 

V||  >  0,  where  Te0  is  the  thermal  electron  temperature  at  1  R©. 
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With  the  kinetic  effects  of  suprathermal  electrons  being  considered,  the  total  elec¬ 
tron  distributions  can  no  longer  be  in  the  form  of  a  drifting  Maxwellian.  The  distri¬ 
butions  are  now  formed  by  the  combination  of  a  drifting  Maxwellian  for  the  thermal 
electrons,  and  the  results  from  the  kinetic  calculations  for  the  suprathermal  electrons. 
As  in  the  basic  version  of  the  model,  the  various  parameters  that  constitute  the  drift¬ 
ing  Maxwellian  distributions  for  the  thermal  electrons  -  as  well  as  the  ambipolar 
electric  field  -  are  determined  in  the  fluid  calculations. 

The  density  and  drift  velocity  of  the  thermal  electrons  are  determined  from  the 
quasi-neutrality  and  current-free  constraints  that  characterize  the  solar  wind,  with 
input  of  densities  and  flow  velocities  of  the  other  particle  components  (proton,  helium, 
and  suprathermal  electron)  obtained  from  the  kinetic  part  of  the  model.  In  particular, 
the  current-free  condition  is: 

neue  +  nsiis  =  npUp  +  2naua,  (III.2.21) 

where  the  subscripts  p,  and  a  stand  for  the  protons  and  the  alpha  particles  respec¬ 
tively.  The  density  and  drift  velocity  profiles  of  the  thermal  electrons  (and  in  fact,  of 
every  plasma  component)  should  also  satisfy  the  continuity  equation: 

I(t)=°-  <iil2-22> 

However,  note  that  both  Eq.  (III.2.21)  and  (III.2.22)  involve  the  number  flux  of 
the  particle  species.  Because  the  profiles  of  the  particle  components  obtained  from 
the  kinetic  part  of  the  model  have  already  satisfied  Eq.  (III. 2. 22),  the  imposition  of 
the  current-free  constraint  throughout  the  whole  simulation  range  renders  thermal 
electron  density  and  velocity  profiles  that  would  also  satisfy  Eq.  (III.2.22). 

The  thermal  electron  temperature,  together  with  the  ambipolar  electric  field,  are 
obtained  from  the  moment  equations  of  the  next  two  higher  orders,  i.e.  the  momen¬ 
tum  and  energy  transfer  equations  for  the  whole  electron  population: 
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where  e  is  the  elementary  charge,  4>g  is  the  gravitational  potential,  is  the  potential 
for  the  field-aligned  ambipolar  electric  field,  Qw>s  is  the  energy  flux  of  the  suprather- 
mal  electron  population,  obtained  from  the  kinetic  calculations,  and  SMe,s/8t  and 
8EetS/8t  represent  the  rates  of  momentum  and  energy  transfers  to  the  thermal  and 
suprathermal  electrons  due  to  Coulomb  collisions.  Unlike  the  standard  moment 
method  [S chunk,  1977],  where  the  collisional  transfer  values  are  completely  based 
on  assumed  distribution  functions,  such  as  a  Maxwellian,  the  rates  of  momentum  and 
energy  transfers  in  Eq.  (III. 2. 23)  and  (III. 2. 24)  are  determined  based  on  the  results 
of  the  kinetic  components  of  our  model.  By  calculating  the  momentum  and  energy 
transferred  to  the  ions  from  the  electrons  in  our  kinetic  simulations,  we  find  these  col¬ 
lisional  transfer  values  with  momentum  and  energy  conservation,  which  results  from 
the  anti-symmetries  of  the  Coulomb  collisional  operator  in  Eq.  (III. 2. 12).  However,  we 
should  note  that  the  suprathermal  electrons  are  not  included  in  the  collisional  back¬ 
ground  for  the  ions.  Nor  are  the  ions  included  in  the  collisional  background  for  the 
suprathermal  electrons.  These  collisional  values,  therefore,  do  not  take  into  account 
the  transfer  between  the  ions  and  the  suprathermal  electrons.  Fichtner  et  al.  [1996] 
studied  the  collisional  transfers  based  on  a  number  of  distribution  functions.  Their 
results  indicated  that  the  transfer  values  for  some  non-Maxwellian  distributions  could 
be  significantly  different  from  those  for  a  Maxwellian.  However,  they  also  noted  that 
for  distributions  whose  core  regions  do  not  deviate  significantly  from  a  Maxwellian, 
their  collisional  transfer  values  are  similar  to  the  Maxwellian  case.  As  we  shall  show, 


the  suprathermal  electrons  in  our  calculations  do  not  significantly  affect  the  core  of 
the  total  electron  distributions.  Therefore,  we  believe  that  the  transfer  values  due  to 
collisions  between  the  suprathermal  electrons  and  the  ions  are  negligible. 

We  note  that  the  formulation  of  the  fluid  calculations  is  simpler  when  the  kinetic 
effects  of  suprathermal  electrons  are  not  considered.  In  that  case,  the  total  electron 
distributions  consist  only  of  the  thermal  electrons,  and,  by  assumption,  are  in  the 
form  of  a  drifting  Maxwellian.  One  may  then  drop  all  the  terms  with  subscript  s 
in  Eq.  (III.2.21),  (III.2.23)  and  (III. 2. 24)  to  recover  the  fluid  equations.  However, 
whether  suprathermal  electron  effects  are  considered  or  not,  when  the  results  of  both 
the  fluid  and  kinetic  parts  of  our  model  converge,  not  only  the  collisional  opera¬ 
tors  in  the  kinetic  calculations  provide  a  self-consistent  background  for  the  Coulomb 
interactions  of  the  ions,  but  the  ambipolar  electric  field  is  also  consistent  with  the 
particle  distributions  themselves.  Convergence  of  the  results  thus  renders  physically 
meaningful  solutions  for  our  solar  wind  model. 

III.2.3  Comparison  of  Solar  Wind  Driving  Mechanisms 
Strong  Wave-Driven  Case 

As  discussed,  we  can  incorporate  the  suprathermal  electron  effects  into  our  solar 
wind  model  by  including  the  optional  kinetic  component  described  in  Section  III.2.2. 
Because  interactions  between  the  electrons  and  the  waves  are  not  considered,  the 
kinetic  suprathermal  electron  effects  in  our  calculations  are  essentially  the  velocity 
filtration  effect  proposed  by  Scudder  and  Olbert  [1979].  Thus,  by  comparing  solutions 
generated  with  and  without  the  optional  suprathermal  electron  component  being 
included,  we  can  identify  the  contribution  by  the  velocity  filtration  effect  in  driving 
the  solar  wind. 

We  have  generated  two  such  solar  wind  solutions  for  the  range  between  1  R© 
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and  1  AU  with  identical  parameters  and  boundary  conditions.  The  interplanetary 
magnetic  field  that  we  use  is  based  on  a  dipole  -f  quadrupole  +  current  sheet  model 
[Banaszkiewicz  et  al.,  1998]: 

2  4.5  1 


B{r)  =  2.21 


+ 


Gauss.  (III.2.25) 


[(r/R0)3  (r/R0)5  1.538(r/R0  +  1.538)2J 

Boundary  and  initial  conditions  are  imposed  at  the  lower  end  of  the  calculations:  a 
thermal  electron  temperature  of  100  eV  is  assumed.  Initial  velocity  distributions  for 
the  protons  and  helium  ions  are  the  upper  halves  of  Maxwellian  distributions,  whose 
temperatures  are  also  100  eV.  The  proton  and  helium  ion  densities  are  respectively 
4.2  x  107  and  4.0  x  106  cm-3.  The  parameter  rj0  that  characterizes  the  strength  of 
the  available  wave  power  is  0.016.  As  we  shall  discuss  below,  such  a  value  for  rj0 
corresponds  to  the  case  of  a  strong  wave-driven  solar  wind. 

Comparing  the  two  solar  wind  solutions,  we  find  that  the  most  noticeably  different 
results  are  the  shape  of  the  total  electron  distributions.  Although  in  both  cases  we 
begin  with  a  total  electron  distribution  in  the  form  of  a  Maxwellian  at  1  R0,  we 
find  that  the  distributions  become  different  downstream  in  the  flow.  Figure  III. 2.1 
shows  the  total  electron  distributions  at  two  different  locations.  At  0.1  AU,  in  the 
solution  where  global  kinetic  effects  are  considered  for  the  suprathermal  electrons,  a 
significant  tail  in  the  outward  portion  of  the  total  electron  distribution  has  already 
been  developed.  Downstream  in  the  solar  wind,  even  though  the  temperature  of  the 
thermal  electrons  is  reduced,  the  tail  remains  as  strong  (middle  panel  of  Fig.III.2.1). 
Thus,  relative  to  the  thermal  core,  the  tail  portion  of  the  electron  distribution  is  more 
enhanced  downstream.  The  formation  and  enhancement  of  this  non-thermal  feature 
verify  Scudder  and  Olbert’s  [1979]  idea  of  the  velocity  filtration  effect. 

The  deviation  of  the  total  electron  distribution  from  the  Maxwellian  approxima¬ 
tion  can  be  seen  more  clearly  in  the  bottom  panel  of  Fig.III.2.1,  where  the  reduced 
distributions  in  the  parallel  direction  for  the  two  solutions  are  shown.  It  is  evident 
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Figure  III-2.1:  Top  and  middle  panels:  Contour  plots  for  the  total  electron  distributions 
at  0.1  AU  (top)  and  0.4  AU  (middle)  for  the  case  that  includes  suprathermal  electrons. 
Bottom  panel:  Reduced  electron  distributions  in  the  parallel  direction  at  0.4  AU.  Solid  line 
is  the  distribution  for  the  case  where  global  kinetic  effects  of  the  suprathermal  electrons  are 
included.  Dashed  line  is  the  result  with  a  Maxwellian  approximation  for  the  entire  electron 
population. 
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Figure  III.2.2:  Solid  (dotted)  line:  Parallel  (perpendicular)  temperature  of  the  total  elec¬ 
tron  population  for  the  case  where  global  kinetic  effects  of  the  suprathermal  electrons  are 
included.  Dashed  line:  Electron  temperature  for  the  Maxwellian  approximation.  The  thick 
line  segment  and  the  cross  (x)  at  the  right  of  the  plot  correspond  to  the  core  electron 
temperatures  observed  by  the  Helios  [Pilipp  et  al.,  1990]  and  Ulysses  [Scime  et  al.,  1994] 
spacecraft  respectively. 

that  while  the  suprathermal  electrons  have  a  very  significant  contribution  in  the  tail 
portion  of  the  distribution,  they  do  not  really  affect  the  core  region.  In  fact,  by 
comparing  the  overall  electron  temperature  in  the  two  solutions,  we  find  that  the 
suprathermal  electrons  do  not  affect  the  results  by  a  significant  amount  (Fig.  III.2.2). 
Therefore,  our  test-particle  approximation  for  the  suprathermal  electrons  regarding 
their  Coulomb  collisions  seems  reasonable. 

The  formation  of  the  suprathermal  tail  in  the  electron  distributions  leads  to  an 
outward  total  electron  heat  flux  throughout  the  solar  wind,  as  opposed  to  a  zero 
heat  flux  with  a  Maxwellian  approximation.  As  shown  in  Fig.  III. 2.3,  the  profile  of 
this  outward  electron  heat  flux  essentially  follows  a  power  law  from  beyond  10  R©  (~ 
0.05  AU).  Although  the  heat  flux  profile  has  a  power-law  index  of  —2.06,  which  seems 
to  agree  with  the  estimates  based  on  a  number  of  fast  solar  wind  events  observed 
by  the  Helios  spacecraft  [Pilipp  et  al .,  1990],  its  value  at  1  AU,  like  the  electron 
temperature,  is  about  an  order  of  magnitude  too  low  compared  with  the  observations 
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Figure  III.2.3:  Profile  of  the  total  electron  heat  flux  from  1  R®  to  1  AU  in  the  case  where 
the  kinetic  effects  of  the  suprathermal  electrons  are  included.  The  meaning  of  the  thick  line 
segment  and  the  cross  (x)  is  similar  to  that  in  Fig.  III. 2.2. 

(see  Fig.  III.2.2  and  III.2.3).  However,  it  has  been  shown  that  kinetic  wave-particle 
interactions  involving  the  solar  wind  electrons  can  increase  the  temperature  and  the 
heat  flux  of  the  species  remarkably  [  Tam  and  Chang ,  2001].  In  the  present  study,  we 
exclude  the  electron  heating  so  as  to  present  a  clearer  picture  of  the  impact  by  the 
velocity  filtration  effect. 

The  suprathermal  electron  population  influences  the  ion  species  mainly  through 
the  ambipolar  electric  field.  The  ambipolar  electric  potential  profiles  for  the  two  solar 
wind  solutions  are  shown  in  Fig.  III. 2. 4.  With  a  Maxwellian  approximation  for  the 
entire  electron  population,  the  overall  electric  potential  drop  is  about  700  V  from  1  R® 
to  1  AU.  The  potential  difference  increases  by  about  70  V  when  kinetic  suprathermal 
electron  effects  are  taken  into  account.  Such  an  increase  suggests  that  the  velocity 
filtration  effect  can  accelerate  the  ions  in  the  solar  wind. 

To  evaluate  the  significance  of  the  velocity  filtration  effect  as  a  solar  wind  driving 
mechanism,  we  shall  compare  its  contribution  with  those  from  other  physical  mecha- 
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Figure  III.2.4:  Ambipolar  electric  potential  profiles.  Solid  line:  kinetic  suprathermal  elec¬ 
tron  effects  included;  dashed  line:  with  a  Maxwellian  approximation  for  the  entire  electron 
distribution. 

nisms  in  terms  of  energy  input  to  the  ion  populations.  Besides  the  kinetic  effect  that 
associates  with  the  suprathermal  electrons,  there  are  several  physical  processes  affect¬ 
ing  the  ion  energy  in  the  solar  wind.  First,  ion  cyclotron  resonant  heating  provides 
a  direct  source  of  energy  input  for  both  the  protons  and  helium  ions.  In  addition, 
Coulomb  collisions  allow  energy  exchange  between  different  species  in  the  solar  wind. 
The  energies  of  the  ions  are  affected  by  these  two  physical  processes  not  only  directly, 
but  also  indirectly  through  the  influence  of  these  mechanisms  on  the  ambipolar  elec¬ 
tric  field.  As  the  ions  travel  downstream  along  the  solar  wind  flow,  their  kinetic 
energy  increases  due  to  the  sunward  gradient  of  the  ambipolar  electric  potential  (see 
Fig.  III. 2.4),  but  at  the  same  time,  some  of  the  kinetic  energy  is  converted  into  gravi¬ 
tational  potential  energy.  Whether  an  ion  can  reach  large  radial  distances  in  the  solar 
wind  depends  whether  the  kinetic  energy  it  gains  is  high  enough  to  overcome  the  grav¬ 
itational  potential.  Therefore,  we  evaluate  the  significance  of  a  physical  mechanism 
in  driving  the  solar  wind  by  comparing  its  contribution  to  the  energy  input  for  the 
ions  with  the  contributions  from  all  the  physical  processes  combined.  Because  we  are 
interested  in  comparing  only  the  processes  that  may  drive  the  solar  wind,  we  shall  ex- 
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Figure  III.2.5:  £  for  protons  and  helium  ions.  Solid  lines:  kinetic  suprathermal  electron 
effects  included;  dashed  lines:  with  a  Maxwellian  approximation  for  the  entire  electron 
distributions. 

elude  the  gravitational  force  in  our  energy  consideration.  After  all,  the  gravitational 
potential  profile  is  invariant  under  all  solar  wind  conditions,  and  is  independent  of 
the  presence  of  other  physical  processes.  Hence,  we  define  the  following  quantity  for 
each  ion  species  for  our  purpose  of  energy  comparison: 

£  =  average  kinetic  energy  +  gravitation  potential  energy  per  ion. 

For  a  given  ion  species  in  our  model,  its  increase  in  £  along  the  solar  wind  flow  is 
due  to  a  combination  of  energy  input  from  the  ambipolar  electric  field,  wave-particle 
interactions,  and  Coulomb  collisions  with  other  species.  In  particular,  the  ambipolar 
electric  field  also  reflects  the  influence  due  to  wave-particle  interactions,  and  Coulomb 
collisions,  and,  in  the  case  where  the  kinetic  suprathermal  electron  effects  are  included, 
the  velocity  filtration  effect  as  well.  The  profiles  of  £  for  the  protons  and  helium  ions 
in  our  two  solar  wind  solutions  are  shown  in  Fig.  III. 2.5.  It  appears  from  the  figure 
that  the  difference  in  £  is  very  small  between  the  two  solutions.  Indeed,  we  find 
that  with  the  kinetic  effects  of  the  suprathermal  electrons  taken  into  account,  the 
increase  in  £  for  the  protons  from  1  R0  to  1  AU  is  only  1.5%  larger  than  in  the  case 
of  a  Maxwellian  approximation.  For  the  helium  ions,  the  corresponding  difference  is 
2.4%. 
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Figure  III.2.6:  Ion  temperature  profiles.  Solid  lines:  kinetic  suprathermal  electron  effects 
included;  dashed  lines:  with  a  Maxwellian  approximation  for  the  entire  electron  distribu¬ 
tions. 

Because  the  quantity  S  characterizes  the  average  kinetic  energy,  it  is  related  to 
both  the  temperature  and  average  outflow  speed  of  the  species.  The  temperature 
profiles  of  the  ions  in  the  two  solar  wind  solutions  are  compared  in  Fig.  III. 2.6.  One 
can  see  from  the  figure  that  the  influence  on  the  proton  and  helium  ion  temperatures 
due  to  the  velocity  filtration  effect  is  negligible. 

The  influence  by  the  kinetic  suprathermal  electron  effects  on  the  ion  outflow  speeds 
can  be  seen  in  Fig.  III.2.7.  The  difference  due  to  the  velocity  filtration  effect  is 
minimal,  considering  that  the  kinetic  effect  only  lead  to  a  1.5%  increase  in  the  proton 
speed  and  a  2.3%  increase  in  the  helium  ion  speed  at  1  AU.  Note  that  even  under  a 
Maxwellian  approximation  for  the  overall  electron  distributions,  the  solar  wind  speed 
at  1  AU  is  as  high  as  650  km/s  in  the  solution.  Thus,  the  available  wave  power 
alone  in  these  calculations  seems  to  be  able  to  drive  the  solar  wind  velocity  to  the 
high-speed  range.  Because  of  that,  these  solutions  correspond  to  the  case  of  a  strong 
wave-driven  solar  wind. 

As  the  ions  in  the  two  solar  wind  solutions  have  similar  energies,  outflow  speeds, 
and  temperatures,  their  velocity  distributions  are  also  similar.  Fig.  III. 2.8  shows  the 
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Figure  III. 2.7:  Profiles  of  the  proton  and  helium  ion  outflow  speeds.  Solid  lines:  kinetic 
suprathermal  electron  effects  included;  dashed  lines:  with  a  Maxwellian  approximation  for 
the  entire  electron  distributions. 

cuts  of  the  normalized  proton  velocity  distributions  at  the  vj_  =  0  plane  of  these 
strong  wave-driven  solutions.  We  can  see  that  the  distributions  of  the  two  solutions 
do  not  show  a  significant  difference  (open  circles  and  crosses  in  Fig.  III.2.8).  For  these 
solutions,  there  are  no  clear  indications  of  a  secondary  peak.  We  have  generated  a 
different  strong  wave-driven  solution  with  kinetic  suprathermal  electron  effects  con¬ 
sidered.  The  parameters  of  this  solution  are  identical  to  the  cases  discussed  above, 
except  that  the  ratio  of  sunward  to  anti-sunward  propagating  wave  power,  instead  of 
Eq.  (III.2.5),  is  replaced  by  the  following  expression: 

e_/e+  =  1  -  e~s2  .  (III.2.26) 

As  shown  in  Fig.  III.2.9,  e_/e+  is  larger  by  about  an  order  of  magnitude  with 
Eq.  (III.2.26)  replacing  Eq.  (III. 2. 5).  The  increase  in  the  sunward  propagating  wave 
power  does  not  seem  to  affect  the  core  region  of  the  proton  velocity  distributions 
by  much,  as  indicated  in  Fig.  III. 2.8.  However,  it  has  a  stronger  effect  on  the  tail 
of  the  distributions.  With  the  increase  e_/e+  ratio,  indications  of  the  formation  of 
a  secondary  peak  at  the  tail  are  much  stronger  (filled  circles  in  Fig.  III.2.8).  This 
result  supports  the  association  between  double-peaked  proton  velocity  distributions 
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Figure  III.2.8:  Cuts  of  normalized  proton  velocity  distributions  at  the  =  0  plane  at 
0.2  AU.  The  crosses  (x)  corresponds  to  the  solution  based  on  a  Maxwellian  approximation 
for  the  total  electron  population.  Open  circles  (o)  are  for  the  solution  that  includes  kinetic 
suprathermal  electron  effects.  The  ratio  e_/e+  is  given  by  Eq.  (III. 2. 5)  for  both  of  these 
cases.  Filled  circles  (•)  are  for  e_/e+  given  by  Eq.  (III. 2.26),  with  suprathermal  electron 
effects  considered.  For  clarity,  the  distributions  are  normalized  differently  so  that  their  peak 
values  are  exactly  one  order  of  magnitude  apart.  The  straight  solid  lines  at  the  tail  of  the 
distributions  are  parallel  and  also  exactly  one  order  of  magnitude  apart. 


Figure  III.2.9:  The  ratio  of  sunward  to  anti-sunward  propagating  wave  power.  Open  circles: 
based  on  Eq.  (III.2.5);  filled  circles:  based  on  Eq.  (III.2.26). 
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and  sunward  propagating  waves  that  was  discussed  earlier. 

Overall,  it  is  clear  from  our  comparison  of  the  ion  energies,  outflow  speeds,  tem¬ 
peratures,  and  distributions  that  when  strong  wave-particle  interactions  are  present, 
the  velocity  filtration  effect  plays  an  insignificant  role  in  the  driving  and  acceleration 
of  the  solar  wind.  Recently,  Vocks  and  Marsch  [2001,  2002]  have  introduced  a  semi- 
kinetic  model  for  the  solar  wind  to  describe  the  effects  of  strong  ion  cyclotron  resonant 
heating.  Like  this  study,  their  model  describes  the  wave-particle  interactions  in  the 
framework  of  quasilinear  theory.  However,  the  treatment  of  the  particles  is  different 
in  the  two  models.  Their  model  is  based  on  reduced  ion  distributions,  massless  elec¬ 
tron  fluid,  and  an  assumed  electron  temperature  profile.  With  the  approximations 
on  the  electrons,  the  model  does  not  take  into  account  the  kinetic  effects  associated 
with  the  suprathermal  electrons.  However,  as  we  have  shown  that  the  velocity  filtra¬ 
tion  effect  associated  with  the  suprathermal  electrons  is  negligible  when  strong  ion 
resonant  heating  is  present,  Vocks  and  Marsch’s  assumptions  regarding  the. electrons 
seem  to  be  reasonable  approximations. 

Weak  Wave-Driven  Case 

Kinetic  suprathermal  electron  effects  were  proposed  to  be  a  solar  wind  acceleration 
mechanism  under  the  assumption  that  no  wave-particle  interactions  were  present 
[Olbert,  1982].  Therefore,  it  is  worthwhile  for  us  to  investigate  the  velocity  filtration 
effect  in  the  scenario  of  a  weak  wave-driven  solar  wind.  In  this  scenario,  we  assume 
that  the  available  wave  power  for  ion  resonances  is  not  sufficient  to  drive  the  outflow 
to  the  velocity  range  of  the  fast  solar  wind.  We  are  interested  in  finding  out  whether 
the  velocity  filtration  effect  can  make  a  considerable  contribution  in  driving  the  solar 
wind  in  this  situation. 

We  compare  two  solar  wind  solutions  for  a  weak  wave-driven  case,  one  with  the 
kinetic  effects  associated  with  the  suprathermal  electrons  taken  into  account,  and  the 
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Figure  III. 2. 10:  Profiles  of  the  proton  and  helium  ion  outflow  speeds  for  the  weak  wave- 
driven  solar  wind  case.  Solid  lines:  kinetic  suprathermal  electron  effects  included;  dashed 
lines:  with  a  Maxwellian  approximation  for  the  entire  electron  distributions. 

other  with  a  Maxwellian  approximation  for  the  entire  electron  distribution.  These 
calculations  essentially  share  the  same  boundary  conditions  as  those  in  the  previous 
section,  but  with  two  exceptions.  First,  rjo  is  reduced  to  0.0075  to  characterize  the 
weak  wave-driven  scenario.  Second,  as  indicated  in  Fig.  III. 2.5,  the  profiles  for  S  have 
plateaued  long  before  the  distance  reaches  1  AU.  Therefore,  for  the  present  purpose, 
it  is  sufficient  for  us  to  set  the  upper  boundary  of  the  calculations  to  0.1  AU  (~ 
21.5  R®). 

With  the  reduced  level  of  available  wave  power,  and  without  the  suprathermal 
electron  effects,  we  find  that  both  the  proton  and  helium  ion  speeds  are  only  about 
390  km/s  at  0.1  AU  (Fig.  III. 2. 10).  The  inclusion  of  kinetic  suprathermal  electron 
effects  increases  these  speeds  only  by  2.5%  for  the  protons  and  2.8%  for  the  helium 
ions. 

With  the  reduced  wave  power,  the  ion  parallel  temperatures,  as  shown  in 
Fig.  III. 2. 11,  are  smaller  by  about  half  compared  with  the  corresponding  values  in 
the  strong  wave-driven  solar  wind  case  (Fig.  III.2.6).  However,  as  in  the  previous 
case,  the  velocity  filtration  effect  does  not  seem  to  affect  the  ion  temperature  in  the 
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Figure  III. 2. 11:  Ion  temperature  profiles  for  the  weak  wave-driven  solar  wind  case.  Solid 
lines:  kinetic  suprathermal  electron  effects  included;  dashed  lines:  with  a  Maxwellian  ap¬ 
proximation  for  the  entire  electron  distributions. 

present  scenario. 

As  for  the  energy  quantity  £,  comparison  between  Fig.  III. 2. 5  and  III.2.12  indicates 
that  the  values  for  both  ion  species  reduce  by  about  one-thirds  due  to  the  weaker  wave 
power.  By  comparing  the  two  solutions  in  Fig.  III. 2. 12,  one  can  see  that  with  the 
reduced  wave  power,  the  velocity  filtration  effect  is  still  only  responsible  for  a  small 
portion  of  the  total  energy  input  to  the  ions.  For  both  ion  species,  inclusion  of  the 


Heliocentric  Distance  (Re) 


Figure  III.2.12:  £  for  protons  and  helium  ions  in  the  weak  wave-driven  solar  wind  scenario. 
Solid  lines:  kinetic  suprathermal  electron  effects  included;  dashed  lines:  with  a  Maxwellian 
approximation  for  the  entire  electron  distributions. 
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kinetic  effect  leads  to  an  increase  of  less  than  2%  in  the  total  energy  input  from  1  R® 
to  0.1  AU.  Therefore,  even  in  a  weak  wave-driven  scenario,  the  velocity  filtration 
effect  is  negligible  in  terms  of  driving  and  accelerating  the  solar  wind. 

III. 2.4  Discussion 

Our  results  have  shown  that  ion  cyclotron  resonant  heating  is  a  much  more  efficient 
solar  wind  driving  and  acceleration  mechanism  compared  with  the  velocity  filtra¬ 
tion  effect  that  associates  with  the  suprathermal  electrons.  Nevertheless,  kinetic 
suprathermal  electron  effects,  in  general,  may  not  always  be  negligible.  In  situations 
where  there  are  physical  processes  which  affect  the  suprathermal  electron  population 
more  strongly  than  Coulomb  collisions,  kinetic  suprathermal  electron  effects  may 
contribute  to  the  solar  wind  acceleration  considerably.  In  fact,  the  strength  of  these 
kinetic  effects  depends  on  the  electron  distribution.  We  note  that  the  suprathermal 
electrons  considered  in  this  study  come  from  the  tail  portion  of  a  Maxwellian  dis¬ 
tribution,  whose  temperature  is  100  eV.  Based  on  our  experience,  a  higher  electron 
temperature,  due  to  the  kinetic  effects  of  the  suprathermal  electrons,  would  translate 
into  a  larger  total  electron  heat  flux.  That  would  lead  to  a  larger  ambipolar  electric 
field,  and  consequently,  higher  ion  outflow  speeds.  Such  a  dependence  between  the 
electron  temperature  and  the  solar  wind  flow  speed  has  been  demonstrated  by  Mak- 
simovic  et  al.  [1997].  But  perhaps  more  importantly,  an  electron  distribution  that 
deviates  from  the  Maxwellian  by  having  an  enhanced  suprathermal  tail  may  affect 
the  solar  wind  speed  even  more  strongly.  Maksimovic  et  al.  [1997]  and  Pierrard  et  al. 
[1999]  studied  the  effects  of  electron  Kappa  distributions  on  the  solar  wind  from  a  few 
R®  to  1  AU.  They  showed  that  in  cases  where  the  electron  distributions  had  a  very 
enhanced  suprathermal  tail,  the  solar  wind  velocities  could  approach  the  observed 
high-speed  range.  Although  no  physical  processes  were  related  to  the  possible  for- 
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mation  of  such  highly  non-thermal  electron  distributions,  these  authors  pointed  out 
the  association  between  a  Kappa-like  distribution  near  the  corona  and  the  observed 
electron  distributions  at  1  AU. 

A  physical  mechanism  that  can  significantly  enhance  the  tail  portion  of  the  solar 
wind  electron  distributions,  and  increase  the  electron  temperature  and  heat  flux  is 
wave-particle  interactions.  Chashei  and  Fahr  [2000]  suggested  that  electron  heating 
due  to,  for  example,  preexisting  whistler  waves  in  the  solar  wind,  is  required  to  pro¬ 
duce  the  electron  heat  flux  observed  by  the  Ulysses  spacecraft.  More  recently,  one  of 
our  studies  [Tam  and  Chang ,  2001]  investigated  the  effects  of  electron  cyclotron  reso¬ 
nant  heating  in  the  solar  wind.  We  showed  that  the  electron  heating  by  the  right-hand 
polarized  waves  near  the  corona  could  raise  the  electron  heat  flux  and  temperature 
by  an  order  of  magnitude,  and  lead  to  a  much  larger  ambipolar  electric  field  due  to 
the  associated  kinetic  effects.  The  enhanced  electric  field  resulted  in  a  16%  increase 
in  the  proton  outflow  speed  at  1  AU  for  the  strong  wave-driven  case.  Thus,  even 
though  ion  resonant  heating  appears  to  be  the  dominant  solar  wind  driving  mech¬ 
anism,  the  kinetic  effects  associated  with  the  suprathermal  electrons  can  sometimes 
make  a  considerable  contribution  to  the  acceleration  of  the  solar  wind. 

III. 2. 5  Conclusions 

In  this  study,  we  have  compared  the  effects  of  two  solar  wind  acceleration  mechanisms: 
wave-particle  interactions  and  the  kinetic  effects  of  the  suprathermal  electrons.  Be¬ 
cause  the  suprathermal  electron  population  in  our  calculations  corresponds  to  the 
tail  portion  of  a  Maxwellian  distribution,  it  has  been  implicitly  assumed  that  these 
suprathermal  electrons  have  not  undergone  any  heating,  such  as  that  involving  elec¬ 
tron  cyclotron  resonances.  In  addition,  besides  Coulomb  collisions,  we  do  not  consider 
any  physical  process  which  may  strongly  affect  the  shape  of  the  electron  distributions. 
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Therefore,  the  kinetic  suprathermal  electron  effects  that  we  consider  are  essentially 
the  velocity  filtration  effect,  which  arises  due  to  Coulomb  collisions  and  the  global 
kinetic  nature  of  the  solar  wind  flow.  Our  results  have  shown  that  in  the  presence 
of  ion  cyclotron  resonant  heating,  the  contribution  by  the  velocity  filtration  effect  in 
the  driving  and  acceleration  of  the  solar  wind  is  relatively  insignificant.  However, 
when  there  are  other  physical  processes,  such  as  electron  cyclotron  resonances,  that 
strongly  affect  the  suprathermal  electrons,  kinetic  suprathermal  electron  effects  may 
no  longer  be  negligible. 
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